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Abstract 

The decay of a false vacuum of unbroken B—L symmetry is an intriguing and 
testable mechanism to generate the initial conditions of the hot early universe. If 
B—L is broken at the grand unification scale, the false vacuum phase yields hy- 
brid inflation, ending in tachyonic preheating. The dynamics of the B—L break- 
ing Higgs field and thermal processes produce an abundance of heavy neutrinos 
whose decays generate entropy, baryon asymmetry and gravitino dark matter. 
We study the phase transition for the full supersymmetric Abelian Higgs model. 
For the subsequent reheating process we give a detailed time-resolved description 
of all particle abundances. The competition of cosmic expansion and entropy 
production leads to an intermediate period of constant 'reheating' temperature, 
during which baryon asymmetry and dark matter are produced. Consistency of 
hybrid inflation, leptogenesis and gravitino dark matter implies relations between 
neutrino parameters and superparticle masses. In particular, for a gluino mass 
of 1 TeV, we find a lower bound on the gravitino mass of 10 GeV. 
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1 Introduction 



Neutrino masses, baryogenesis, dark matter and the acoustic peaks in the power spec- 
trum of the cosmic microwave background (CMB) radiation require an extension of the 
Standard Model of particle physics. The supersymmetric standard model with right- 
handed neutrinos and spontaneously broken B—L, the difference of baryon and lepton 
number, provides a minimal framework which can account for all these phenomena [l]. 
B—L breaking at the grand unification (GUT) scale leads to an elegant explanation 
of the small neutrino masses via the seesaw mechanism and explains baryogenesis via 
leptogenesis [2] . The lightest supersymmetric particle is an excellent candidate for dark 
matter pHs] and the spontaneous breaking of B—L requires an extended scalar sector, 
which automatically yields hybrid inflation [6|[7] , explaining the inhomogeneities of the 
CMB. 

Recently, we have suggested that the decay of a false vacuum of unbroken B—L 
symmetry generates the initial conditions of the hot early universe: nonthermal and 
thermal processes produce an abundance of heavy neutrinos whose decays generate 
primordial entropy, baryon asymmetry via leptogenesis and gravitino dark matter from 
scatterings in the thermal bath (sIIq]. In this context, tachyonic preheating after hybrid 



inflation 10 sets the stage for a matter dominated phase whose evolution is described by 
Boltzmann equations, finally resulting in a radiation dominated phase. It is remarkable 
that the initial conditions of this radiation dominated phase are not free parameters but 
are determined by the parameters of a Lagrangian, which in principle can be measured 
by particle physics experiments and astrophysical observations. 



Our work is closely related to previous studies of thermal leptogenesis 11,12 and 



nonthermal leptogenesis via infiaton decay 13 16 , where the infiaton lifetime deter- 
mines the reheating temperature. In supersymmetric models with global B—L symme- 
try the scalar superpartner Ni of the lightest heavy Majorana neutrino A^^i can play the 



role of the infiaton in chaotic 17 ,18 or hybrid 19 ,20 inflation models. One of the main 
motivations for nonthermal leptogenesis has been that the 'gravitino problem' for heavy 



unstable gravitinos 21-25 can be avoided by means of a low reheating temperature. In 
the following we shall assume that the gravitino is the lightest superparticle. Gravitino 
dark matter can then be thermally produced at a reheating temperature compatible 



with leptogenesis 26 
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The present work is an extension of Ref. |9]. We discuss in detail the effect of all 
supersymmetric degrees of freedom on the reheating process and restrict the param- 
eters of the Lagrangian such that they are compatible with hybrid inflation and the 
production of cosmic strings during spontaneous symmetry breaking. This implies in 
particular that B—L is broken at the GUT scale. The consistency of hybrid inflation, 
leptogenesis and gravitino dark matter entails an interesting connection between the 
lightest neutrino mass nii and the gravitino mass m^. As we shall see, the final results 
for baryon asymmetry and dark matter are rather insensitive to the effects of super- 
particles and details of the reheating process. Due to the restrictions on the parameter 
space compared to Ref. j9] the lower bound on the gravitino mass increases to about 
10 GeV. 

The paper is organized as follows. In Section [2] we briefly recall field content and su- 
perpotential of our model, in particular the Froggatt-Nielsen flavour structure on which 
our analysis is based. We then discuss the time-dependent masses of all particles dur- 
ing the spontaneous breaking of B—L symmetry in the supersymmetric Abelian Higgs 
model, the restrictions of hybrid inflation and cosmic strings on the parameters, and the 
particle abundances produced during tachyonic preheating. Section [3] deals with the 
time evolution after preheating and the required set of Boltzmann equations for all par- 
ticles and superparticles. The detailed description of the reheating process is given in 
Section |4] with emphasis on the various contributions to the abundance of A^i neutrinos, 
the lightest of the heavy Majorana neutrinos, whose decays eventually generate entropy 
and baryon asymmetry. Particularly interesting is the emerging plateau of a reheating 
temperature which determines the final gravitino abundance. In Section [5] a system- 
atic scan of the parameter space is carried out, and relations between neutrino and 
superparticle masses are determined. Three appendices deal with important technical 
aspects: The full supersymmetric Lagrangian for an Abelian gauge theory in unitary 
gauge, which is used to describe the time-dependent B—L breaking (Appendix [A| , 
CP violation in all supersymmetric 2 — )■ 2 scattering processes (Appendix [b]) and the 
definition of the reheating temperature (Appendix [C| . 
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2 B—L breaking at the GUT scale 



2.1 Field content and superpotential 

Our study is based on an extension of the minimal supersymmetric standard model 
(MSSM) which offers solutions to a series of problems in particle physics and cosmology. 
Its main features are right-handed neutrinos, a U{1)b-l factor in the gauge group and 
three chiral superfields, needed for B—L breaking and allowing for supersymmetric 
hybrid inflation. In this section, we give a review of this model, presented earlier in 
Ref. [9], thereby focussing on the aspects which are especially relevant for this paper. 

A characteristic feature of the model is that inflation ends in a phase transition which 
breaks the extra U{1) symmetry. During this phase transition the system experiences 
the decay from the false into the true vacuum. At the same time, this phase transition 
is responsible for the production of entropy, matter and dark matter through tachyonic 
preheating and subsequent leptogenesis. Finally, it yields masses for the right-handed 
neutrinos, thereby setting the stage for the seesaw mechanism, which can explain the 
observed light neutrino masses. The superpotential is given by 

W = ^^ (^l-L - 2 3,82) + -^hMn^S, + ht^^n^H^ + W^mssm , (1) 

where Si and 5*2 are the chiral superfields containing the Higgs field responsible for 
breaking B—L, $ contains the infiaton, i.e. the scalar field driving inflation, and n1 
denote the superfields containing the charge conjugates of the right-handed neutrinos. 
In the following, we will refer to the components of 5*1, S2 and $ as the symmetry 
breaking sector, whereas the components of n1 form the neutrino sector. Vb-l is the 
scale at which B—L is broken. The B—L charges are qs = qs2 = —QSi = 2, g$ = 0, 
and gn, = —1- h and A denote coupling constants, and PFmssm represents the MSSM 
superpotential, 

W^MSSM = h'l^lOaO.H^ + hf^bllO.H^ . (2) 

For convenience, all superfields have been arranged in SU (5) multiplets, 10 = (g, u'^, e^) 
and 5* = (rf^, /), and i, j = 1, 2, 3 are flavour indices. We assume that the colour triplet 
partners of the electroweak Higgs doublets if„ and have been projected out. The 
vacuum expection values f„ = {H^) and Vd = (Hd) break the electroweak symmetry. In 
the following we will assume large tan /3 = Vu/vd, implying Vd Vu — wew = ^/f^ + v^. 
For notational convenience, we will refer to Hu as H from now on. 
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Table 1: Froggatt-Nielsen flavour charge assignments. 



In addition to these chiral superfields, the model also contains a vector supermul- 
tiplet V ensuring invariance under local B—L transformations and the gravity super- 
multiplet consisting of the graviton G and the gravitino G. 



2.2 Froggatt-Nielsen flavour model 

The flavour structure of the model is parametrized by a Froggatt-Nielsen flavour model 



based on a global f/(l)FN group, following Refs. 27,28]. According to this model, the 
couplings in the superpotential can be estimated up to 0{1) factors as powers of a 
common hierarchy parameter rj, with the exponent given by the sum of the flavour 
charges Qi of the fields involved in the respective operators. Setting the charges of all 
Higgs fields to zero, this implies 

hij ~ r^^'+'^j ^ ~ r7<3* . (3) 

The numerical value of the parameter ~ 1 / -\/300 is deduced from the quark and lepton 
mass hierarchies. This remarkably simple flavour model can reproduce the experimental 
data on Standard Model masses and mixings, while at the same time it remains flexible 
enough to incorporate the phenomena beyond the Standard Model mentioned above. 
Further details on the predictive power of this model can be found in Ref. |29], where 
we recently performed a Monte-Carlo study to examine the impact of the 0{1) factors. 

In the following, we will restrict our analysis to the case of a hierarchical heavy 
(s)neutrino mass spectrum. Mi ^ M2,M3, where M = K^vb-l- Furthermore we 
assume the heavier (s)neutrino masses to be of the same order of magnitude as the 
common mass mg of the particles in the symmetry breaking sector, for definiteness 
we set M2 = M3 = ms- With this, the Froggatt-Nielsen flavour charges are fixed as 
denoted in Tab. [!} Taking the B—L gauge coupling to be g"^ = g^ux — 7r/6, the model 
can now, up to 0{1) factors, be parametrized by the f/(l)FN charges a and d. The 
B—L breaking scale vb-l, the mass of the lightest of the heavy (s)neutrinos Mi, and 
the effective light neutrino mass parameter mi are related to these by 

Ml vb-l 



Here, fn^ = y/m2m3, the geometric mean of the two hght neutrino mass eigenvalues m2 
and m3, characterizes the hght neutrino mass scale, which, with the charge assignments 
above, can be fixed to 3 x 10~^ eV. To obtain this result, we exploited the seesaw 
formula nii, = —mi:>M~^nr[) with m/5 = h'^VEw- Furthermore, it can be shown that rhi 
is bounded from below by the lightest neutrino mass mi (30| . 

In the following, we will study the model in terms of the more physical quantities 
vb^l and Mi instead of the ?7(1)fn charges. To partly account for the 0{1) uncertain- 
ties in the neutrino mass matrices, we will additionally vary rhi. Apart from this, we 
ignore any further uncertainties of the model and simply set the 0{1) prefactors to one. 
Furthermore, when considering the production of dark matter in form of gravitinos, cf. 



Section 3.2, the gravitino (m^) and gluino {rrig) masses will be additional parameters. 



2.3 Spontaneous symmetry breaking 

Before the spontaneous breaking of B—L, supersymmetry is broken by the vacuum 
energy density po = I-^^b-l? which drives infiation. During this time, the dynamics 
of the system is governed by the slowly rolling scalar component (p of infiaton 
multiplet $. The scalar components of the Higgs superfields 5*1 2 are stabilized at zero. 
The right-handed sneutrinos and the scalar MSSM particles obtain their masses due 
to supergravity contributions. As the field value of the infiaton decreases, so do the 
effective masses in the Higgs sector, until a tachyonic direction develops in the effective 
scalar potential. The subsequent phase transition can best be treated in unitary gauge, 
in which the physical degrees of freedom are manifest. In particular, performing a 
super-gauge transformation relates the Higgs superfields 5*1,2 and the vector superfield 
V to the respective fields S' and Z in unitary gauge, 

S,,2 = ^S' exp{±tT), V = Z+^^{T-T*). (5) 

Note that the chiral superfield T playing the role of the gauge transformation parameter 
is chosen such that Si and 5*2 are mapped to the same chiral superfield S'. This refiects 
the fact that one chiral superfield is 'eaten' by the vector superfield in order to render 
it massive. The supermultiplet S' contains two real scalar degrees of freedom, s' = 
;^(a'+2r), where r remains massive throughout the phase transition and a' is the actual 
symmetry-breaking Higgs field. It acquires a vacuum expectation value proportional 
to v{t) = ^(cr'^(t, x))y^ which approaches vb-l at large times. In the Lagrangian, we 
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account for symmetry breaking by making the replacement a' — \/2v{t) + a, where a 
denotes the fluctuations around the homogeneous Higgs background. 

The fermionic component s of the supermultiplet S' pairs up with the fermionic 
component of the inflaton supermultiplet $ to form a Dirac fermion the higgsino, 
which becomes massive during the phase transition. Due to supersymmetry, the corre- 
sponding scalar fields (a, r and inflaton 0) end up having the same mass as the higgsino 
in the supersymmetric true vacuum. Likewise, the gauge supermultiplet Z (gauge bo- 
son A, real scalar C, Dirac gaugino A) and the (s)neutrinos Ni {Ni) acquire masses. 
Note that the choice of unitary gauge, cf. Eq. (|5]), forbids us to use the Wess-Zumino 
gauge, so Z denotes a full massive gauge multiplet with four scalar and four fermionic 
degrees of freedom. The capital N refers to the physical Majorana particle N = {n, n)'^ 
built from the two Weyl spinors contained in the superfields n*^ and n. N denotes the 
complex scalar superpartner of the left-chiral fermion n. For an overview of the particle 
spectrum, see Fig. [Tj 

At the end of the phase transition, supersymmetry is restored. An explicit calcula- 
tion of the Lagrangian describing this phase transition is given in Appendix lAl We can 



read off the mass eigenvalues during the phase transition, cf. Eqs. (106) to (110): 

ml = lx{3v\t) - , ml = ^X{vI_l + v\t)) , 

ml = Xv^{t), ml = Xv^{t), 
ml = %gW{t) , 
Ml = {h^fv\t) . 

Here we have ignored corrections which arise due to thermal effects and due to super- 
symmetry breaking before the end of inflation in some hidden sector, leading to a mass 
for the gravitino. 



2.4 Hybrid inflation and cosmic strings 

The spontaneous breaking of B—L discussed in the previous section marks the end of 
a stage of hybrid inflation, which is governed by the flrst term in the superpotential in 
Eq. ([T|. As the symmetry breaking proceeds very rapidly and abruptly, it represents 
what is often referred to as a 'waterfall' phase transition. It is accompanied by the 
production of local topological defects in the form of cosmic strings as well as the 
nonadiabatic production of particles coupled to the Higgs fleld, a process commonly 
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Figure 1: Nomenclature, production and decay processes after B~L breaking. The Higgs field a and 
particles coupled to it are produced during tachyonic preheating, as marked by the red boxes. The 
gauge degrees of freedom then decay nearly instantaneously (black, dashed arrows), whereas the decay 
and production of the other degrees of freedom can be described by Boltzmann equations (blue, solid 
arrows). The numbers in parentheses denote the respective internal degrees of freedom. 



known as tachyonic preheating 10 . In the following, we shall first discuss cosmic 



strings and then tachyonic preheating. 

Due to the nontrivial topology of its vacuum manifold, the Abelian Higgs model 
underlying the B—L phase transition gives rise to solitonic field configurations. These 
are called cosmic strings (for a review, cf. e.g. j^)- Their energy per unit length is 

(7) 



32 



the 



/x = 27ri?(/3K_^, 

with /3 = A/(8 and B{l3) = 2.4 [ln(2//3)]~^ for (3 < lO'^. According to Ref. 
characteristic length separating two strings formed during tachyonic preheating is 

e = (-A^;B-L0c)-'/^ (8) 
Here 0c is the velocity of the radial component of the inflaton field, = (/9/-\/2e*^, at 



9 



the onset of the phase transition, which can be determined from the scalar potential by 
exploiting the equation of motion for ip in the slow-roll approximation, 3iJ0 = —V'{(p). 



According to Refs. 33 , 34 , in the region of parameter space we are interested in. 



the slope of the scalar potential is determined by the Coleman- Weinberg one-loop 



corrections, cf. e.g. Ref. 35 for the explicit formulas resulting from the superpotential 
Eq. g. With this, the energy density stored in strings just after the end of the phase 
transition can be calculated as 

Pstring = ^ ■ (9) 

Inserting Eqs. ([T]) and ^ and the expression for the Coleman- Weinberg one-loop poten- 
tial from Ref. [35], we find that the fraction of energy stored in cosmic strings directly 
after the phase transition increases strongly with the coupling parameter A. This is due 
to the higher energy density per cosmic string as well as the shorter average distance 
between two strings. For instance, for vb-l = 5 x 10^^ GeV and A = 10"^, we find 
{HE,y^ ~ 400 and Pstring/po ^ 60 %. For A = 10-^ this is reduced to {H^Y^ ~ 40 and 
Pstring /Po =^0.2%. 

As the universe evolves, the cosmic strings intercommute, forming closed loops which 
are separated from the infinite strings. These oscillate, loosing energy into gravitational 
waves as well as into the Higgs and gauge degrees of freedom until they eventually decay, 
cf . Refs. (36|[37] . After a relaxation time, which is roughly given by tstring ~ ^ [36] , there 
is only 0{1) cosmic string per Hubble volume left and the energy density stored in the 
cosmic strings scales as Pstring oc H'^Mp. These relic cosmic strings can in principal 
be observed today, e.g. via string induced gravitational lensing effects in the CMB. 
The nonobservation of these effects implies an upper bound on the energy per unit 



length 38 41 . In the following, we will work with 



Gp<5xl0-^ (10) 

where G = Mp"^ is Newton's constant with Mp = 1.22 x 10^^ GeV denoting the Planck 
mass. Inserting this into Eq. ([T]) puts an upper bound on Vb^l, which weakly depends 
on A, 

^B-L<1.8xl0-Mln^) Mp. (11) 



In Ref. 35 , the authors discuss hybrid inflation and cosmic string production in a 
setup very similar to oursj^ Taking into account current experimental bounds inferred 
^Cf. also the analyses in Refs. [34 ' 



42 



43 
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from the spectrum of fluctuations in the CMB |44 and from the nonobservation of 



cosmic strings 38 , they find viable inflation for 

3 X 10^^ GeV < vb-l < 7 x 10^^ GeV 



(12) 

10-^ < < 10-1 . 



This significantly constrains the allowed parameter space. With the scale of B—L 
breaking basically fixed, vb-l ~ 5 x 10^^ GeV, Eq. Q implies a = and a factor of 
proportionality of about 5. This is still consistent with the Froggatt-Nielsen model, 
since three 0{1) factors enter in the calculation of vb-l- The bounds on A restrict 
the second free t/(l)FN charge, 1.4 < d < 2.6, cf. Eq. ([3]), and therefore Mi. In the 
following, we will consider the restricted parameter space 

Vb-l = 5 X 10^^ GeV, 

10^ GeV < Ml < 3 X 10^^ GeV , (13) 
10"^ eV < mi < 1 eV. 

Here, the variation of fhi accounts for the uncertainties of the Froggatt-Nielsen model. 



The chosen range easily covers the expected values for rhi in this setup, cf. Ref. 29 
for a recent analysis. 

The production and decay of cosmic strings can in principle have a large influence 
on the state of the universe just after the phase transition. However, as we will argue in 
the following, for our purposes it is not necessary to treat these processes in detail, as 



long as we restrict ourselves to the parameter space in Eq. (13). According to Eq. ([8j), 
it is possible to have as much as (9(10^) cosmic strings per Hubble volume for large 
values of the coupling parameter A. For the maximal value of the coupling constant, 
A = 10"^, roughly half of the total energy density just after the phase transition is 
stored in cosmic strings. However, since in this case the relaxation time of the cosmic 
strings, tstring ~ — O {10^^) , is much smaller than a Hubble time, the major 
component of this energy has been converted back into Higgs and gauge degrees of 
freedom before the processes which we describe by means of Boltzmann equations, cf. 



Section 3.2 , become relevant. Since the exact mechanism of energy loss of cosmic strings 
is not yet fully understood, we will in the following omit the effects from cosmic strings, 
keeping in mind that at the very most, they will convert about half of the initial energy 
density of the Higgs bosons into particles of the Higgs and gauge multiplets. Typically 
the effects from cosmic strings are much less important, e.g. for A < lO^'' their relative 
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energy contribution is at the level of at most (9(1%). Due to supersymmetry, the 
extra higgsinos produced will decay into the same supermultiplet as the Higgs bosons 
would have, thus inducing no significant change in the following discussion. The extra 
gauge particles will decay predominantly into radiation, which is quickly diluted at 
this early stage of the matter dominated phase governed by the nonrelativistic Higgs 
bosons. Thus it can be expected that our setup is insensitive to a modification of the 
contribution from cosmic strings to the initial conditions of the following reheating 
phase. We also confirmed this in a numerical study. Considering the case of extremal 
string production, we shifted half of the energy initially stored in the Higgs bosons at 
the end of preheating into the gauge degrees of freedom and calculated the resulting 
entropy, baryon asymmetry and gravitino dark matter. We find no deviations from the 
results presented in Section |4] above the percent level. 



2.5 Tachyonic preheating 

Let us now consider the production of particles coupled to the Higgs field. As the value 
of the infiaton field decreases, the scalar potential develops a tachyonic instability 
in the direction of the Higgs field. Quantum fiuctuations of the Higgs field cr^ with 
wave number \k\ < |mo-| begin to grow exponentially, while its average value remains 
zero. The strong population of the long wavelength Higgs modes leads to a large 
abundance of nonrelativistic Higgs bosons. Other particles coupled to the Higgs field 



are nonperturbatively produced due to the rapid change of their effective masses 45 



The mode equations for the gauge, Higgs, infiaton, and neutrino supermultiplets 
are governed by the time-dependent masses proportional to v{t) given by Eq. 
According to Ref. 45 , this leads to particle production, with the energy and number 
densities for bosons and fermions after tachyonic preheating given 



Pb/po ~ 2 X 10 ^gsX /(a, 1.3) , nij(a) ~ 1 x 10 ^^.m|/(a, 1.3)/a , 

p^/po ^ 1.5 X 10-3(^,A/(a,0.8), n^(a) ~ 3.6 x 10-^(7,m|/(a, 0.8)/a , (14) 

with /(a, 7) = a/q;^ + 7^ — 7 and a = mx/ms, where mx denotes the mass of the 
respective particle in the true vacuum; Qs counts the spin degrees of freedom of the 
respective particle. Just as the Higgs bosons themselves, these particles are produced 
with very low momentum, i.e. nonrelativistically. 



^Note that particle production can be significantly enhanced by quantum effects 46 , which require 
further investigations. 
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A deviation from this mechanism is found for the imaginary component r of the 
complex field s', since fixing the gauge to unitary gauge yields a constant contribution 
to its mass. Neglecting the expansion of the universe, the mode equation for r reads 

d^n + {k'' + ml)Tk = . (15) 

We can absorb the constant mass contribution in the momentum k. In the language of 
Ref. ^5|, this is equivalent to a shift in the 'asymptotic in frequency' uj-{k). To excite 
a given mode more energy is necessary, the production is thus less efficient. 

With the ingredients discussed so far, the stage is now set for the emergence of 
the hot early universe. In other words, B—L breaking after hybrid inflation provides 
the initial conditions for the successful generation of entropy, matter and dark matter. 
Before we demonstrate this numerically, we will first introduce the necessary tools, i.e. 
decay rates and supersymmetric Boltzmann equations, in Section [3] 



3 Time evolution after preheating 
3.1 Decay rates and branching ratios 

In this section, we will discuss the evolution of the particle abundances from the initial 
conditions set by B—L symmetry breaking to the radiation dominated era. A schematic 
overview of all relevant processes is given in Fig. [TJ 

During tachyonic preheating, most of the vacuum energy is converted into Higgs 
bosons. At the same time, particles coupled to the Higgs bosons, i.e. particles of 
the gauge, Higgs, infiaton and neutrino supermultiplets are produced (cf. red boxes 
in Fig. [T]), with the resulting abundances given by Eq. (14). Among these particles. 



the members of the gauge supermultiplet have by far the shortest lifetime. Due to 
their large couplings they decay basically instantaneously into (s)neutrinos and MSSM 
particles (cf. dashed, black arrows in Fig. [T]). This sets the initial conditions for the 
following phase of reheating, which we will describe with Boltzmann equations (cf. 
solid, blue arrows in Fig. [T]). 

Due to our choice of a hierarchical (s)neutrino mass spectrum, the decay of parti- 
cles from the symmetry breaking sector into the two heavier (s) neutrino generations is 
kinematically forbidden. These particles can hence only decay into particles of the A^^i 
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supermultiplet. These (s) neutrinos, just as the neutrinos produced through gauge par- 
ticle decays and thermally produced (s)neutrinos, decay into MSSM particles, thereby 
generating the entropy of the thermal bath as well as a lepton asymmetry. Note that 
these different production mechanisms for the (s)neutrinos yield (s)neutrinos with dif- 
ferent energies, which due to relativistic time-dilatation, decay at different rates. Fi- 
nally, the thermal bath produces a thermal gravitino abundance, which will turn out 
to be in the right ball-park to yield the observed dark matter abundance. 

In the following we list the total and partial vacuum decay rates necessary to quan- 
tify the processes described above. The total vacuum decay rates for the particles of 
the symmetry breaking, gauge and neutrino sectors are 



r".-r2,..* = 3^W)^".(i-4^) . (16) 

with X denoting the superfields of the model carrying B—L charges qx- The relevant 
partial decay rates at leading order are given by 

po _ _ _ 4 ^ po po - (i- aMTX fO 

po = 1 = 1 = i 

A->-0x0x 2 A^i^xi^x 2 C^(t>x(t>x 2 A-^(j)xipx 



1 



(19) 



/ 2 \ 1/2 



with (f)x and ipx denoting the scalar and fermionic components of a superfield X. 
At tree level the pseudoscalar r decays exclusively into fermionic neutrinos, similar 
to its scalar partner a, whose branching ratio into scalar neutrinos is suppressed by 
two powers of the mass ratio Mi/ms- The production of r particles during tachyonic 



preheating, cf. Sec. 2.5, is however negligible compared to the production of a particles. 
We can thus neglect the contribution from the pseudoscalar r in the following. 
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Another important consequence of Eq. (19) is that taking into account that the 



production of the different gauge degrees of freedom during tachyonic preheating is 



proportional to the respective spin degrees of freedom, cf. Eq. (14), the decay products 
of the gauge degrees of freedom consist to equal shares of scalar and fermionic degrees 
of freedom. 

3.2 Super symmetric Boltzmann equations 

In this section, we present the Boltzmann equations describing the evolution of the 
universe after the decay of the gauge degrees of freedom, as depicted by the blue, solid 
arrows in Fig. [TJ This analysis is a supersymmetric extension of the study performed 
earlier in Ref. |9], exploiting the techniques explained there in detail. In general, the 
evolution of the phase space density fxit,p) of a particle species X is determined by 

^fx{t,p) = Y,Y.^x{Xah..^i3.:), (20) 

ab.. ij.. 

with C denoting the Liouville operator and the Cx denoting the collision operators of 
all relevant processes involving the particle X: 

^fx{t,p) =jjx{t,p) , 

Cx{Xah..^ij..)=-^y] [ dU{X\a,b, .:){2n)^6^''\Pout - P^n) (21) 

^3xEx ^ J 

X [/./,.. |-M(U.. ^ Xab..)\' - fxfah..\M{Xah.. ^ , 



where 'Ylidof denotes the sum over all internal degrees of freedom of the initial and final 
states and the momentum space element cill is given by 

cl?p 

cin(X|a,6, ..) = S{X,a,b, ..^ij, ..) dpadpb..dpidpj.. , dp = j^^^^^j^ ■ (22) 

S{X, a, b, ..; 2, j, ..) is a statistical factor to prevent double counting of identical particles. 
The quantum statistical factors due to Bose enhancement and Pauli blocking have been 
omitted, since typically they only yield minor corrections |47) . 

In the following, we will often work with integrated Boltzmann equations, which 



are obtained by integrating Eq. (20) over gx d^px/i^n)^ . In a Friedmann-Lemaitre- 



Robert son- Walker Universe, the resulting equation can be simplified to 

aH^Nx = TxNx , (23) 
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with a denoting the scale factor, H the Hubble rate, Tx the effective production rate 
of X particles, and 



a{t) 



9x 



(24) 



GeV J (27r)3 

the comoving number density, i.e. the number of X particles in a volume (a/GeV)'^. 



Performing a rescaling of a in Eq. (23) leaves the physical number density nx invariant. 
For convenience, we will thus set apH = 1 at the end of preheating in the following. 
Another useful quantity next to the number density nx is the energy density px, 



px = 9x 



:Ex{p)fx{t,p). 



(25) 



Evolution of the gravitational background 

The time-dependence of the scale factor a{t) is governed by the Friedmann equation. 
For a flat universe and a constant equation of state u = p/p between some reference 
time to and time t, the Friedmann equation yields 



a(t) = aito) 



1+2(1 + ^) 



Stt 
3M|' 



1/2 

Ptot(to) ) {t - to) 



(26) 



After preheating, the universe is dominated by nonrelativistic Higgs bosons, i.e. u = 0. 
After the end of the reheating process, the universe is radiation dominated, u = 1/3. In 
the intermediate region, the equation of state changes continuously. We approximate 
this by implementing a piecewise constant effective equation of state with coefficients Ui 
in the intervals (ti,ti^i] with a^jj < a(ti) < a(tj+i) < a^^jj. The Ui are determined 
iteratively by requiring self-consistency of the Friedmann equation, 

3(1+Wi) 

(27) 



Ptot{ti 



a(t 



i+l) 



Ptot{ti+l) V (^(ti) 

In our numerical calculations, we approximate the total energy density by its two 
dominant components, the energy density of the Higgs bosons and the energy density 
of the neutrinos produced in Higgs, higgsino and infiaton decays, ptot 



Pa + P%,: 



for which we will obtain analytical expressions below, cf. Eqs. (35) and (37). In the 



following we will calculate the Hubble rate H = a/a using Eq. (26). 
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Massive degrees of freedom 

The Boltzmann equations describing the massive degrees of freedom introduced above 
are 

tf. = -C,{u ^ iViiVi) - C,{a ^ iViiVi) , (28) 

tf<i, = -C^{<p ^ NiNi) , (29) 

CU = -C4tlj^N*Ni), (30) 

CJn,,, = -C^,,3(iV2,3^MSSM), (31) 

= -C^^,,3(iV2,3^MSSM), (32) 

CfN, = 2 CnM NiNi) + CN,{iJ N;Ni) + Cn,{Ni^ MSSM), (33) 

+ C^^(iVi ^ MSSM). (34) 
a, and ip, the particles of the symmetry breaking sector, are produced via tachyonic 



preheating only, hence their initial number densities are given by Eq. (14) and their 



initial phase space distributions are peaked at low momenta, and thus can be taken to 



be proportional to S{p). The collision operators on the right-hand side of Eqs. (28) to 



(30) describe the decay of these particles. The resulting ordinary differential equations 



are solved by 



fx{t,p) = — iV^(tpH)7^ exp[-r^ - ^ph)] , X = a, 0, , (35) 
9x {apY 

with tpH denoting the time at the end of preheating. We fix the origin of the time axis 
by setting tpH = 0. Also the abundances of all heavy (s)neutrinos obtain contributions 
from tachyonic preheating. The corresponding phase space distribution functions are 



of the same form as fx in Eq. (35). 



The collision operators for the lightest (s) neutrinos are more involved. Just as in 
Ref . [9j , they can be treated best by separating the phase space density into the contri- 
butions due to thermal (th) and nonthermal (nt) (s)neutrinos. Introducing Sx{Eq] t, t'), 
the energy of a particle X at time t which was produced with energy time t', 



ox[Eo; t ,t) = Eq 



a{t) 



a{t) 




(36) 



17 



we find for tlie comoving number densities of nontliermally produced (s)neutrinos: 



iV-(t) = iV^(t) + iV^^(t) + iV|(t) 



(37) 



dt' a^(t') 'js,xit') exp 



dt"- 



Ni 



+ iVPH(tpH)e- 



Sx{ms/2;t',t")\ 

r /"* 

%it-tPu)^ NG^tc) exp - / dt' 

L Jta 

witli X = A^i, iVi and 



(38) 



Here A^™(tpH) denotes tlie initial X abundance from nonperturbative particle produc- 
tion during tachyonic preheating, whereas Nxitc) refers to the initial X abundance 
from the decay of the gauge degrees of freedom. Note that this notation is valid through- 
out this paper: The lower indices on number densities, decay rates, etc. indicate the 
respective particle, whereas the upper index refers to the origin of this particle. The 
time to denotes the lifetime of the gauge particles after preheating, tc = tpH + '^/'^% cf. 
Eq. (16), and corresponds to the value ac of the scale factor, ac = a (tc)- Furthermore, 



also the (s)neutrinos of the second and third generation are produced in the decays of 
gauge particles. The corresponding comoving number densities of these (s)neutrino 



species are of the same form as in Eq. (37) 



Inserting Eq. (35) and Eq. (19) into Eq. (38) yields the same time-dependence, 
'-ys,x oc exp(— r° (t — tpn)), for both neutrinos and sneutrinos. Eq. (37) hence implies a 



constant ratio between neutrinos and sneutrinos produced via decays of particles from 
the symmetry breaking sector throughout the reheating phase. For instance. 



ATS" 



4.5 X 10"^ 



for Ml = lO^^GeV. 



(39) 



The precise value and the dependence on Mi arises due to the initial conditions set by 



tachyonic preheating, cf. Eq. (14), and the branching ratios denoted in Eq. (19). For 
increasing Mi, we find a weak increase of A^^ /^ivi- 

Unlike the two heavier (s)neutrino generations, (s)neutrinos of the first generation 
are also produced thermally from the bath. Assuming kinetic equilibrium, their co- 
moving number densities are determined by the integrated Boltzmann equation 



da 



X 



-r^(A^^^ - A^; 



X 



X = ATi, N, , 



(40) 
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with N'^ denoting the comoving number density in thermal equihbrium and r^_^j^- is 
the vacuum decay width weighted with the average inverse time dilatation factor, 



- - r° [ Ml p (ao\ 



In Eq. (40) we are interested in the decay width of the thermally produced neutrinos, 



r*^^. In this case Eq. (|42|) can be evaluated to T% = Ki{Mi/T) / K2{Mi/T), where 
Kn denotes the modified Bessel function of the second kind of order n. Note however 



that Eq. (|42|) is not restricted to this case but also allows the calculation of, for example. 



the decay width of the neutrinos produced by the decay of the Higgs bosons, Ff 



MSSM degrees of freedom 

The Boltzmann equations governing the lepton number asymmetry and the abundance 
of MSSM particles in the thermal bath are 

tfL = C, + C~,-Ci-Ci,, (43) 
tfR = r{C, + C-, + Ci + C-,.), (44) 

with Cf^i denoting the collision operators responsible for the production, decay and 
scattering of (anti)(s)leptons and r describing the number of radiation quanta produced 
in the respective processes. 

A subtle but important point concerning the Boltzmann equation for the lep- 
ton asymmetry is the correct treatment of 2 — > 2 scattering processes with heavy 
(s)neutrinos in the intermediate state. The collision operator for (s)neutrino decay 
takes care of the on-shell contributions to these processes, so we need to add the off- 
shell contributions. The CP-conserving part is negligible compared to the on-shell 
contribution, so we shall concentrate on the CP-violating part. This can be obtained 
by calculating the CP-violating contribution of the full 2 — )■ 2 scattering process and 
then subtracting the on-shell CP- violating contribution (reduced collision operator). 
By exploiting unitarity and CPT invariance, we prove in Appendix [B] that the CP- 
violating contribution of the full 2 — i- 2 scattering process vanishes up to corrections 
of 0{{h^Y), so that we can replace the CP-violating off-shell contribution by the neg- 
ative of the CP-violating on-shell contribution. With this, the integrated Boltzmann 



equation up to (9(ej, {h'^Y) obtained from Eq. (43) reads 

d 



(^H^Nl = rf iVf + T'tN'^ - TwNl (45) 
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with the washout rate tw and the effective (non)thermal production rates for the lepton 



asymmetry r^'"^* given by 



i-w = — J- 



f?-(ivf)"'E E (r^^AT^H + r^ATf + riivl) , (46) 

In Eq. (45) we have introduced A'^£* and as the nonthermal and thermal contri- 
butions to the total lepton asymmetry = + Nf^, respectively. The decay rate 
of the thermally produced (s)neutrinos, F^^, as well as the decay rates F^^, F^, and 
F^ for nonthermally produced (s)neutrinos, are given by Eq. (42). Note that Eq. (46) 
relates decay rates F and effective production rates f . The latter describe the relative 
increase of the respective particle species due to a given production process and can 
directly be compared with the Hubble rate H in order to determine the efficiency of 
the respective process, parametrizes the CP asymmetry in the Ni and Ni decays, 



which, in the Froggatt-Nielsen model, can be estimated as 48 , 49 



6i < 0.1 . (47) 



In the following, we will set to its maximal value, thus obtaining an upper bound for 
the produced lepton asymmetry. 



Analogously, this time neglecting terms of C(ej), Eq. (44) yields the integrated 
Boltzmann equation for the relativistic degrees of freedom of the MSSM, 

aH^N^ = tj^N^' + t'iN]^ , (48) 
with f denoting the effective rates of (non)thermal radiation production, 

« X=N„N, 



(49) 



Here denotes the effective increase of radiation quanta in the thermal bath by adding 
a particle X stemming from the production mechanism x with energy e^, 



x-TT^- (50) 
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Another important quantity in this context is the total radiation production rate Tr. 
It counts the radiation quanta produced per unit time and is obtained by dividing the 



right-hand side of the Boltzmann equation for radiation, Eq. (48), by A^^^;, 



R 



aH 



~ ill 



(51) 



da " ur i\r 

Here in the last expression, denotes the number density of the dominant source for 
radiation production at a given time. 



Solving Eq. (48) finally yields the temperature T of the thermal bath, 

1/3 



T 



TT 



Nr. 



(52) 



.9*,nC{5) a^ 

with g^:^n counting the effective relativistic degrees of freedom contributing to the num- 
ber density of the thermal bath hr, in the MSSM g^^n = 427/2. 

Gravitinos 

Gravitinos are predominantlj]^ produced through supersymmetric QCD scattering pro- 
cesses in the thermal bath. The corresponding integrated Boltzmann equation is 



aH^Nr 
da ^ 



(53) 



In QCD, up to leading order in the strong gauge coupling ^f^, the effective production 
rate Tq is given by 50 



fG(T) 



{a/GeVy 



N, 



1 + 



a 



3m% 



'■Ml 



rpb 



In 



T2 



m2(T) 



+ 0.8846 



(54) 



with the energy dependent thermal gluino mass, gluon mass and strong coupling con- 
stant 

^'^^^ m^(/io) , mg{T) = ^g,{T)T, 

(55) 



m~g{T) 



3 

8^2 



1 + :r^^s(/^o)ln 



MT) 



-1/2 



with the typical energy scale during reheating estimated as the average energy per 
relativistic particle in the thermal bath, /^(T) ~ ~ 3T. The gravitino mass rriQ and 
the gluino mass at the electroweak scale rrig = mg{^Q) remain as free parameters. 



■^Note that due to the high temperatures reached in this setup, we do not expect a significant 
contribution from nonthermal gravitino production 
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4 The reheating process 



Combining our initial conditions with the Boltzmann equations derived in the previous 
section poses an initial-value problem. Its solution allows us to quantitatively describe 
the generation of entropy, matter and dark matter due to the production and decay of 
heavy (s) neutrinos. We have numerically solved this problem for all values of the input 



parameters within the ranges specified in Eq. (13). In this section we first illustrate 
our findings for a representative choice of parameter values. In Section |5] we then turn 
to the investigation of the parameter space. 

In this paper, we take into account all (super)particles involved in the reheating 
process, in particular the gauge degrees of freedom, which were omitted in earlier 
studies (Sjjo]. This allows us to give a realistic, time-resolved description of the reheating 
process. Furthermore, compared to Refs. [8}[9], we consider a higher scale of B—L 
breaking, vb-l = 5 x 10^^ GeV, which is compatible with hybrid inflation and cosmic 



strings, cf. Section 2.4 However, many of the techniques employed when solving the 
Boltzmann equations are very similar to those discussed in detail in Refs. [8|[9]. We 
hence in the following focus on the physical results, referring the reader to these earlier 
works for more information on the technical aspects. 

4.1 Particle masses and couplings 

Let us study the evolution of the universe after inflation for 

Ml = 5.4 X 10^° GeV, mi = 4.0 x IQ-^ eV , mg = 100GeV, m-g = ITeV . (56) 



As we will see later in Section 5^, requiring successful leptogenesis as well as the right 
gravitino abundance to explain dark matter typically forces Mi to be close to 10^^ GeV. 
Here, we adjust its explicit numerical value such that, given the values for rhi and mg, 
the gravitino abundance comes out right in order to account for dark matter. The 
choice for fhi represents the best-guess estimate in the context of the Froggatt-Nielsen 
flavour model employed in this work, a result we recently obtained in a Monte-Carlo 



study, cf. Ref. 29 . In scenarios of gauge or gravity mediated supersymmetry breaking 
the gravitino often acquires a soft mass of (9(100) GeV, which is why we set mg to 



100 GeV. A gluino mass of 1 TeV is close to the current lower bounds from ATLAS 51 



and CMS |52j. The values in Eq. (56) readily determine several further important 
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model parameters: 



S.OxlO^GeV, (57) 
5.3 X 10"^ 

Here, we have chosen opposite signs for the CP parameters ei and £2,3, so that the sign 
of the total lepton asymmetry always indicates which contribution from the various 
(s) neutrino decays is the dominant one. 

Fig. |2] presents the comoving number and energy densities of all relevant species as 
functions of the scale factor a. In both panels of this figure some of the displayed curves 
subsume a number of closely related species. These combined curves are broken down 
into their respective components in the two panels of Fig. |3] and in the lower panel 
of Fig. |4} The upper panel of Fig. |4] presents the temperature of the thermal bath as 
function of a. In what follows, we will go through the various stages of the evolution 
depicted in Figs. |2||3] and |4] step by step. Subsequent to that we will, based on the plots 
in Fig. [5} discuss the impact of supersymmetry and the particles of the gauge sector on 
our results. 

4.2 Decay of the massive particles 
Initial conditions 

Tachyonic preheating transfers the bulk of the initial vacuum energy into Higgs bosons, 
p<7 (oph) / Po — 1-0, and only small fractions of it into nonrelativistic higgsinos, infla- 
tons, gauge degrees of freedom and (s)neutrinos {Nf^,Nf^). The particles in the 
gauge multiplet decay immediately afterwards around a = ac, giving rise to relativistic 
(s)neutrinos {N^, Np) and an initial abundance of radiation which thermalizes right 
away. Initially, this thermal bath neither exhibits a lepton asymmetry, nor are there any 
gravitinos present in it. The expansion of the universe between the end of preheating 
and the decay of the gauge degrees of freedom is practically negligible, — ctpH = 1- 
Note that technically all plots in Figs. |2} [3] and |4] start at a = a^. 



ms = 1.6 X 10"GeV, 
= 1.9 X lOGeV, 
A = 1.0 X 10-^ 



M2,3 = 1.6 X lO^'GeV, 
r^^3 = 2.1 X 10^°GeV, 
€2,3 = -1-6 X 10-3, 
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Inverse temperature Mi/T 




10" 10' 10^ 10^ 10'* 10^ 10* 10^ 10^ 

Scale factor a 



10 



Inverse temperature M^/T 
10" 10' 



10^ 10' 




Scale factor a 

Figure 2: Comoving number densities (upper panel) and comoving energy densities (lower panel) 

for particles from the symmetry breaking sector (Higgs a + higgsinos ip + inflatons 4>), (non)thermally 
produced (s)neutrinos of the first generation {Nf^+Nl^, A^f'+A^f'), (s)neutrinos of the first generation 
in thermal equilibrium {2N'^'^, for comparison), (s)ncutrinos of the second and third generation {N2.3 + 
^2,3), the MSSM radiation (R), the lepton asymmetry (B—L), and gravitinos (G) as functions of the 
scale factor a. The vertical lines labeled a|^jj, arh and a^jj mark the beginning, the middle and the 



end of the reheating process. The corresponding values for the input parameters are given in Eq. (561 
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Decay of the (s) neutrinos of the second and third generation 

Among all particles present at a = ac the heavy (s)neutrinos of the second and third 



generation have the shortest lifetimes, cf. Eq. (57). Due to time dilatation, the relativis- 
tic (s)neutrinos stemming from the decay of the gauge particles decay slower than the 
nonrelativistic (s)neutrinos produced during preheating. The decay of the (s)neutrinos 
of the second and third generation is consequently responsible for an increase in the 
radiation number and energy densities on two slightly distinct time scales. 

The gauge particles decay in equal shares into neutrinos and sneutrinos, cf. Sec- 
tion |3.1[ Their number densities thus behave in exactly the same way, explaining the 
overlapping curves in Fig. |3| The production of radiation through the decay of these 
iV^3 neutrinos and iV^g sneutrinos is efficient as long as the effective rate of radiation 
production Tji, cf. Eq. (51), exceeds the Hubble rate H. At a ~ 11 it drops below 
the Hubble rate, which roughly coincides with the value of the scale factor at which 
the comoving energy density of radiation reaches its ffist local maximum. The period 
between preheating and this ffist maximum of the radiation energy density can be re- 
garded as the ffist stage of the reheating process. In the following we shall refer to it 
as the stage of A^2,3 reheating. 

Decay of the particles of the symmetry breaking sector 

The production of higgsinos and inflatons during preheating is roughly equally efficient, 
N^p (ciph) /N^ (ciph) — 1-0. Taking into account kinematic constraints resulting from the 



mass spectrum described in Section [272] all particles from the symmetry breaking sector 
exclusively decay into relativistic (s)neutrinos of the ffist generation {Nf,Nf). 

The majority of Higgs bosons, higgsinos and inflatons survives until ts = tpH + l/r^, 

5 



cf. Eq. (16), which corresponds to a scale factor of as — 7.2 x 10 . Roughly up to this 
time the main part of the total energy is stored in these particles. At later times, i.e. 
for a > as, the energy budget is dominated by the energy in radiation]^ Higgs bosons 
that decay earlier than the average lifetime are responsible for the generation of sizeable 
abundances of neutrinos and sneutrinos. The contributions from higgsino and 
inflaton decays to this process are essentially negligible. 



"^Note that in general the value of the scale factor at which the energy in radiation begins to 
dominate is determined by the lifetime of the most long-lived particle. In the case under study the 
Higgs bosons have the longest lifetime, but for other parameter choices this may be the (s)neutrinos 
of the first generation. 
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Inverse temperature Mi/T 
0.05 0.1 0.2 0.5 




Scale factor a 
Inverse temperature Mi/T 




Scale factor a 



Figure 3: Breakdown of the comoving number densities shown in the upper panel of Fig. [2j The 
(s)neutrinos of the second and third generation (iV2^3+A'2,3) (upper panel) split into (s)neutrinos 
that are produced during preheating [N^]^, ^2^) and in the decay of the gauge degrees of freedom 
(^Fsi ^2^3)- all four cases the sum of the contributions from both generations is shown. The 
(s)neutrinos of the first generation (iVf*+A^f*, Nl^+Nf^) (lower panel) split into (s)neutrinos that 
are produced during preheating (A^f ^, -/V™), in the decay of the gauge degrees of freedom (iVp, Ni), 
in the decay of the particles from the symmetry breaking sector {Nf , Nf), and from the thermal bath 
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Production and decay of the nonthermal (s) neutrinos of the first generation 

The decay of the particles from the symmetry breaking sector is the most important 



source for nonthermal (s)neutrinos. According to our discussion in Section 3.2, the ratio 
between the number densities of neutrinos and iVf sneutrinos is fixed to a constant 
value at all times, cf. Eq. (39). For our choice of parameters we find N'? /N§^ ~ 
4.4 X 10~^. Moreover, the large hierarchy between the two decay rates and F^, cf. 



Eq. (57), renders the and number densities unable to exceed the number density 
of the Higgs bosons. From the perspective of the rather long-lived Higgs bosons the 
(s) neutrinos essentially decay right after their production. As long as they are efficiently 
fueled by Higgs decays, the (s)neutrino number densities continue to rise. But once the 
supply of Higgs bosons is on the decline, they die out as well. The overall timescale of 
our scenario is hence controlled by the Higgs lifetime. However, as we will see below, 
the characteristic temperature of the reheating process is by contrast associated with 
the lifetime of the iVf neutrinos. 

Further contributions to the abundances of nonthermal (s)neutrinos come from pre- 
heating as well as the decay of the gauge particles. Just as in the case of the sec- 
ond and third (s)neutrino generation, the nonrelativistic (s)neutrinos produced during 
preheating decay at the fastest rate and the number densities of neutrinos and 
iVf sneutrinos are always the same. 

4.3 Reheating and the temperature of the thermal bath 
Reheating through the decay of neutrinos 

The energy transfer from the nonthermal (s)neutrinos of the first generation to the 
thermal bath represents the actual reheating process. It is primarily driven by the 
decay of the neutrinos which soon have the highest abundance among all (s)neutrino 
species. In analogy to the notion of N2,3 reheating, we may now speak of A^^i reheating. 



This stage of reheating lasts as long as Tr > H, cf. Eq. (51). Let us denote the two 
bounding values of the scale factor at which f = H hy a^jj and a^^. In the case of 
our parameter example we find a^jj ~ 5.3 x 10^ and ~ 9.8 x 10^. Between these 
two values of the scale factor the comoving number density of radiation roughly grows 
like Nji oc a^. Around the comoving energy density of radiation reaches a local 

minimum and around local maximum. Similarly, we observe that the end 

of reheating nearly coincides with the time at which the energy in radiation begins to 
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dominate the total energy budget, a^-^ ~ as- 
Plateau in the evolution of the temperature 

The upper panel of Fig. |4] displays the temperature of the thermal bath T calculated 
according to Eq. (52) as function of the scale factor a. As a key result of our analysis 
we find that during A^^i reheating the temperature stays approximately constant. For 
a between a^j^ and a(^^ it varies by less than an order magnitude. We thus conclude 
that in the first place a^jj and represent the limiting values for a plateau in the 
evolution of the radiation temperature. The origin of this plateau is the continuous 
production of neutrinos during reheating. As long as these neutrinos are produced 
much faster than they decay, their comoving number density grows linearly in time, 
^ Itpu'^^'^ Taking into account that until a ~ the expansion of 

the universe is driven by the energy in the Higgs bosons, i.e. nonrelativistic matter, 
this translates into N^_^ oc a^^^. The number density in turn controls the scaling 
behaviour on the right-hand side of the Boltzmann equation for radiation during A^^i 
reheating, cf. Eq. (48). Using H oc a~'^/^, we find 

<H ^ a < : "^^^^ °^ °^ ""^^^ ' « «^ ^ ~ const. (58) 

Reheating temperature 

The temperature at which the plateau in Fig. |4] is located sets the characteristic tem- 
perature scale of reheating. In addition, it represents the highest temperature that is 
ever reached in the thermal bath as long as one restricts oneself to times at which it 
contains a significant fraction of the total energy budget of the universe, cf lower panel 
of Fig. [2j We define the reheating temperature Trh as the temperature of the thermal 
bath at a = orh, where orh denotes the value of the scale factor when the decay of 
the neutrinos into radiation is about to become efficient. This is the case once the 
Hubble rate H has dropped to the effective decay rate F;^^ , 

(flRH) = H (orh) , Trh = T (orh) . (59) 

This yields a value which is representative for the temperature plateau, cf. Fig. |4j For 
the chosen set of parameters this equation has the following solution, 

aRH^1.9xlO\ if = ~ 3.5 X lO^GeV, Trh ~ 6.1 x 10^ GeV . (60) 

In Figs. [2] and [5] as well as in the upper panel of Fig. |4] the three values of the scale 
factor marking the initial (orh); characteristic intermediate (cirh) and final {a(^^) point 
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Figure 4: (Upper panel) Temperature of the thermal bath T and (lower panel) comoving number 
densities for the nonthermal (-/V£*) and thermal [Nf^) contributions to the total lepton asymmetry 
as well as all (s)neutrino species (iVf* + TVf', Nf" + Nf", 2Nl'^ for comparison and 7V2,3 + ^2,3) as 
functions of the scale factor a. The vertical lines in the upper panel labeled a|^jj, aRH and a^jj mark 
the beginning, the middle and the end of the reheating process. The vertical lines in the lower panel 
respectively mark the changes in the signs of the two components of the lepton asymmetry. 
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of the reheating process are indicated by dashed vertical hnes. A comparison of our 
definition of the reheating temperature with other common approaches can be found 
in Appendix [C} 

Evolution of the temperature away from the plateau 

During N2^3 reheating the temperature first increases up to a maximal value and then 
decreases like a~^/^. The initial rise reflects the production of radiation through decays 
of (s) neutrinos of the second and third generation while the expansion of the universe 
is negligible. The subsequent decrease then follows from the Boltzmann equation for 



radiation, cf. Eq. (48), using the fact that its right-hand side stays almost constant up 



to the end of iV2,3 reheating, 

aG<a< 11: aH^NR (x N^^^^ ^ const. , Nr (x a^/\ T (x a'^'^ . (61) 

Finally, we note that between the two stages of reheating and after the end of reheating 
the temperature drops off like . This is the usual adiabatic behaviour indicating 
that no radiation, i.e. entropy is being produced, 

11 ~ ~ '^RH '^RH ~ ^ • ^ ^ ' ^ const. , T oc . (62) 

4.4 Small departures from thermal equilibrium 

Production and decay of the thermal neutrinos of the first generation 

Unlike the two heavier (s)neutrino flavours the (s)neutrinos of the first generation are 
also produced thermally (A^**^, A^*'^). Thanks to supersymmetry the evolution of the 
A^J^*^ and A^*'^ number densities is governed by exactly the same Boltzmann equation, cf. 



Eq. (40), so that they are identical at all times. As both species inherit their momentum 
distribution from the thermal bath, they are always approximately in kinetic equilib- 
rium]^ Simultaneously, the interplay between decays and inverse decays drives them 
towards thermal equilibrium. Initially, there are no thermal (s)neutrinos present in the 
thermal bath and inverse decays result in a continuous rise of the thermal (s)neutrino 
number densities until a ~ aj^jj. Around this time the temperature drops significantly 
below the mass Mi and the thermal (s)neutrinos become nonrelativistic. The equi- 
librium number density A^^ begins to decrease due to Boltzmann suppression until it 
almost reaches the actual number density of thermal (s)neutrinos. The production of 

^For a more detailed discussion cf. Appendix B of Ref. joj 
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thermal (s)neutrinos can then no longer compete with the expansion of the universe 
and their comoving number densities do not continue to grow. 

This picture, however, soon changes because reheating sets in. As the temperature 
remains almost perfectly constant until a ~ Qrh, the equilibrium number density A?^^ 
is not diminished due to Boltzmann suppression any further up to this time. Instead it 
bends over and starts to increase like the volume, A"^ oc . The number densities of 
the thermal (s)neutrinos subsequently follow this behaviour of the equilibrium number 
density. During the second phase of A^i reheating the temperature slightly decreases 
again, thereby reinforcing the Boltzmann factor in A^^ . Consequently, the equilibrium 
number density stops growing and shortly afterwards starts declining exponentially. 
An instant after it has passed its global maximum, the number densities of the thermal 
(s)neutrinos overshoot the equilibrium number density. Due to their numerical prox- 
imity the two values of the scale factor at which A^^ and A^^^ respectively reach their 
global maxima cannot be distinguished from each other in Figs. |2} Both events occur 
close to a = 6.6 x 10"^. 

Generation of the baryon asymmetry 

The out-of-equilibrium decays of the heavy (s)neutrinos violate L, C, and CP, thereby 
generating a lepton asymmetry in the thermal bath. A first nonthermal asymmetry is 
introduced to the thermal bath during A'2,3 reheating. For ^ a ^ 2.2, the decay of 
the (s)neutrinos stemming from preheating leads to an increase of the absolute value of 
the comoving number density Nf^. In the interval 6.6 ^ a ^ 13 the lepton asymmetry 
is slightly augmented through the decay of the (s) neutrinos which were produced in 
the decay of the gauge particles. The main part of the nonthermal asymmetry is, 
however, generated during A^i reheating, while the scale factor takes values between 
a ~ 2.0 X 10^ and a ~ 1.3 x 10^. At all other times the effective rate at which 
the nonthermal asymmetry is produced is at least half an order of magnitude smaller 
than the Hubble rate. Among all nonthermal (s)neutrinos of the first generation only 
the neutrinos contribute efficiently to the generation of the asymmetry. Their 
decay results in a positive nonthermal asymmetry that gradually overcompensates the 
negative asymmetry produced during A'2,3 reheating. At a ^ 4.6 x 10'^ the entire initial 
asjTiimetry has been erased and A^£* changes its sign. 

Washout processes almost do not have any impact on the evolution of the nonther- 



mal asymmetry. The rate Tw at which these processes occur, cf. Eq. (46), is always 
smaller than the Hubble rate if by a factor of at least (9(10). On top of that, at the 
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time tw is closest to H, which happens around a ~ 4.0 x 10^ when tw/H ~ 0.12, 
the production rate FJ* is constantly larger than Tw by a factor of (9(10), so that the 
effect of washout on the nonthermal asymmetry is indeed always negligible. 

The decays and inverse decays of thermal (s)neutrinos of the first generation are 
responsible for the emergence of a thermal, initially negative asymmetry in the bath. 
As long as the abundance of thermal (s)neutrinos is far away from the one in thermal 
equilibrium, the absolute value of this asymmetry increases rapidly. Around a ~ a^jj 
this is not the case anymore, causing the production of the thermal asymmetry to stall 
for a short moment. At a 6.3 x 10^ the washout rate Tw overcomes the production 
rate f of the thermal asymmetry and its absolute value begins to decline. Note that 
at this time the rates and Tw are smaller than H by roughly a factor 9. Shortly 
afterwards, at a ^ 6.6 x 10^, the number density of thermal (s)neutrinos overshoots the 
equilibrium density which results in the asymmetry being driven even faster towards 
zero. Already at a ~ 2.3 x 10^ the initial thermal asymmetry is completely erased. 
Meanwhile, washout effects recede in importance. From a ^ 6.9 x 10^ onwards, F^'^ 
permanently dominates over tw, which is why, once the thermal asymmetry has turned 
positive, it does not decrease anymore. Instead it freezes out at its maximal value 
around a ^ 4.5 x 10^ which corresponds to the time when the ratio of F^*^ and the 
Hubble rate H drops below 10"^^^. 

The final values of Nf^^ and allow us to infer the present baryon asymmetry 77^ as 
well as its composition in terms of a nonthermal [t]'^) and a thermal (?7^) contribution. 



r)0 0° Af'^t'th 
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(63) 
_ 



Here, Cgph = 8/23 denotes the sphaleron conversion factor, g^^g = 915/4 and 
43/11 stand for the effective numbers of relativistic degrees of freedom in the MSSM 
that enter the entropy density sr of the thermal bath in the high- and low-temperature 
regime, respectively, and = g^/g^^nNn is the comoving number density of photons. 
As final value for the scale factor we use a/ ^ 1.9 x 10^ which is the maximal value 
depicted in the two plots of Fig. [2j In our parameter example we find 

^7B~3.7xlO~^ r/^* ~ 3.7 X 10"^ r^^^ ~ 1.9 x 10"^^ (64) 



Recall that in Section 3.2, we set the CP asymmetry parameter ei to its maximal 



value. In this sense, the resulting values for the baryon asymmetry must be interpreted 
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as upper bounds on the actually produced asymmetry and are thus perfectly compatible 
with the observed value for the baryon asymmetry, rj'^^ ~ 6.2 x 10~^° 44 . We also 



point out that, in fact, the Froggatt-Nielsen model typically predicts values for ei 
that are smaller than the maximal possible value by roughly a factor of (9(10), cf. 
Ref. [29]. Using a generic value for ei according to the Froggatt-Nielsen model rather 
than estimating ei by means of its upper bound, would thus yield an excellent agreement 
between prediction and observation in the context of this parameter example, tjb — rj'^'^. 

Furthermore, we find that in the case under study it is the nonthermal contribution 
rj'^ that lifts the total baryon asymmetry rjB above the observational bound. The 
thermal contribution r/^ is smaller than 77^* by five orders of magnitude. If we discarded 
the entire idea of nonthermally produced (s)neutrinos being the main source of the 
lepton asymmetry and resorted to standard thermal leptogenesis, we would struggle to 
reproduce the observed asymmetry. For the chosen value of mi, standard leptogenesis 
would result in r^^ ~ 10^^'^ which is almost an order of magnitude below the observed 



value, cf. Ref. 12 for details. By contrast, it is still much larger than our result for 
rj^B- This has mainly two reasons. First, in our scenario the decays of the nonthermal 
neutrinos continuously increase the entropy of the thermal bath, cf. Figs. |2] and |4| 
which results in a nonstandard dilution of the thermal asymmetry during and after its 
production. Between, for instance, a ~ 6.3 x 10^, which corresponds to the time when 
the production of the negative asymmetry is reversed and the absolute value of the 
asymmetry starts to decline, and a = aj, the entropy of the thermal bath increases 
by a factor of (9(100). Second, in consequence of the specific reheating mechanism 
at work the generation of the thermal asymmetry is delayed in time, so that it takes 
place at a lower temperature than in the standard case. This implies a correspondingly 
smaller abundance of thermal (s)neutrinos, rendering our thermal mechanism for the 
generation of an asymmetry less efficient. We will resume this comparison of the thermal 
asymmetry r/^ with the expectation from standard leptogenesis in Section 5.2, where 
we will discuss the respective dependence on the neutrino mass parameters rhi and Mi. 

Production of gravitino dark matter 

Inelastic 2 — 2 scattering processes in the supersymmetric thermal plasma, mediated 
predominantly via the strong interaction, are responsible for the production of dark 
matter in the form of gravitinos. As the right-hand side of the gravitino Boltzmann 



equation, cf. Eq. (53), scales like a T , the efficiency of gravitino production in the 
course of reheating is directly controlled by the interplay between the expansion of the 
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universe and the evolution of the temperature. 

During A^2,3 reheating the temperature roughly declines as T oc a"^/^, cf. Eq. (61) 
such that in first approximation 



da ^ 



TqNq oc a^T^ ^ const. , oc oc a'^^^ , A^^ oc a^/^ . 



G 



G 



(65) 



Once the decay of the (s) neutrinos of the second and third generation has ceased, 



the temperature decreases adiabatically, T oc a ^ or equivalently a^T^ oc 



cf. 



Eq. (62 ). The rate of gravitino production Eg then begins to decrease much faster than 
the Hubble rate, in fact, initially even slightly faster than a~^, causing the comoving 
gravitino number density to approach a constant value. The first stage of gravitino 
production is completed around a ~ 28 which corresponds to the time when Eg is half 
an order of magnitude smaller than H. Erom this time onwards. Eg scales like a~^, the 
production term in the Boltzmann equation is negligibly small and is constant. 

The decline of Eg is reversed as soon as the temperature plateau characteristic 
for the phase of Ni reheating is reached such that approximately a^T^ oc a^. While 
Eg <^ H, the gravitino density continues to remain constant and Eg increases al- 
most as fast as a'^. At a ^ 1.9 x 10^ it has nearly caught up again with the Hubble rate, 
i.e. the ratio Eg/-f^ reaches again a value of 10~^/^. This time marks the beginning of 
the second stage of gravitino production. The production term in the Boltzmann equa- 
tion cannot be neglected any longer and, assuming for a moment an exactly constant 
temperature during A^^i reheating, we have 



da ^ 



EgA/'g oc a^T^ oc a^ Eg oc oc a'^^^ 



A^g oc . (66) 



The gravitino density A"g hence begins to grow again, now even faster than during A^2,3 
reheating. This terminates the rise of the rate Eg, turning it into a decline proportional 
to a~^/^. We thus obtain the interesting result that, although the temperature evolves 
differently during A^2,3 and A^i reheating, the rate Eg always runs parallel to the Hubble 
rate during these two stages of the reheating process. 

At the end of A^^i reheating gravitino productions fades away in the same way as 
at the end of A^2,3 reheating. Around a ^ 3.5 x 10^, when Tq/H drops below lO^^''^, 
the gravitino abundance freezes out. The final value of A'g then allows us to calculate 
fl^h^, the present energy density of gravitinos p~ in units of Pc/h'^, 



n„/,2^ rg ^ ^G^'als iVg 
^ Pc/h'^ Pc/h'^ g*,s 



(67) 
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where pc = 1.052 x 10 ^ h"^ GeV cm ^ denotes the critical energy density of the universe, 
h = 0.72 the Hubble rate H in the units H = hx lOOkms"^ Mpc"\ = 410 cm 



the number density of the CMB photons, and g^:^s, 9^ ^^'^ explained below 

Eq. (63). Recall that after fixing mi, mg and rrig we adjusted the heavy neutrino mass. 
Ml = 5.4 X 10^° GeV, such that we would obtain the right abundance of gravitinos to 
account for the observed amount of dark matter Q'^^h'^ ~ 0.11 [44]. By construction, 
we thus now find in our parameter example 

%/i2~0.11. (68) 

In conclusion, we would like to emphasize the intriguing simplicity of this mech- 
anism for the generation of dark matter. Let us in particular focus on the physical 
picture behind the second stage of gravitino production. Initially, at the onset of A'^i 
reheating, the rate Fg is still very small compared to the Hubble rate H. But given 
the constant spacetime density of gravitino production 7g = T^/n^ oc during A''i 
reheating and the rapid growth of the spatial volume due to the expansion, Fg rapidly 
grows sufficiently large to get the production of gravitinos going. During the remaining 
time of A'^i reheating this production can then proceed without further hindrance as 
the universe, although it is expanding, is filled by a thermal bath at a constant temper- 
ature. The continuous production of radiation nullifies the expansion and gravitinos 
are produced as in a static universe. In other words, one key feature of our scenario of 
reheating is that it turns the universe into a chemistry laboratory in which the tem- 
perature is fixed at a certain value so that dark matter can be cooked in it just to the 
right point. 



4.5 Robustness against theory uncertainties 

In the previous part of this section we discussed in detail the emergence of the hot 
thermal universe after infiation. The successful explanation of reheating as well the 
generation of matter and dark matter by means of our scenario did, however, not rely 
on any fortunate coincidence between certain particulars but was a direct consequence 
of the overall setup that we considered. The essential steps in the evolution after 
symmetry breaking were the following. Preheating results in an initial state whose 
energy density is dominated by nonrelativistic Higgs bosons. These decay slowly into 
nonthermal neutrinos of the first generation which in turn decay into radiation, thereby 
reheating the universe, generating a lepton asymmetry and setting the stage for the 
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Inverse temperature Mi/T 




Scale factor a 



Figure 5: Conioving number densities after omitting all massive superparticles (upper panel) and 
in addition the B—L vector boson (lower panel), to be compared with the result of the full anal- 
ysis in Fig. [2] The individual curves show the comoving number densities of the Higgs bosons (a), 
nonthermally and thermally produced neutrinos of the first generation (A^"', N^^), neutrinos from the 
first generation in thermal equilibrium (iV^''), neutrinos of the second and third generation (A'2,3), 
the MSSM radiation (R), the lepton asymmetry [B—L), and gravitinos (G) as functions of the scale 
factor a. The vertical lines labeled o^H' '^RH and a^^i mark the beginning, the middle and the end of 



the reheating process. The corresponding values for the input parameters are given in Eq. ( 56 ) 
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thermal production of gravitinos. At the same time, an additional contribution to 
the lepton asymmetry is generated by thermally produced (s)neutrinos. All further 
details that we took care of are, of course, important for a complete understanding 
of the physical picture, but merely have a small impact on the final outcome of our 
calculation. In particular, as we will illustrate in this subsection, the numerical results 
for the observables of interest, Trh, Vb, and fl^h"^, remain unaffected if one neglects 
the superpartners of all massive particles or if one excludes the gauge particles from the 
analysis, cf. Fig. [5} in which we plot the corresponding comoving number densities of all 
remaining species as functions of the scale factor. This observation renders our scenario 
of reheating robust against uncertainties in the underlying theoretical framework and 
opens up the possibility to connect it to other models of inflation and preheating as 
long as these provide similar initial conditions as spontaneous B—L breaking after 
hybrid inflation. In addition to that, the robustness of our scenario justifies to crudely 
simplify its technical description. If one is interested in the parameter dependence of 
the observables and less in the exact evolution during reheating, one may simply omit 
effects due to the gauge degrees of freedom and supersymmetry as it has been done in 
Refs. fs] and [9]. 

Nonsupersymmetric analysis including the gauge multiplet 

In a first step, in order to assess the impact of supersymmetry on the reheating process 
in the Abelian Higgs model, we neglect the superpartners of all massive particles, i.e. 
the gauge scalar C, the gaugino A, the higgsino ip as well as all heavy sneutrinos Ni. 
Technically, this renders the infiaton (p stable as it can only decay into a pair of A^^i 
sneutrinos. To avoid overclosure of the universe we thus also omit the infiaton. By 
contrast, we keep the full particle spectrum of the MSSM and the gravitino because 
we still wish to account for dark matter by thermally produced gravitinos. All in all, 
these simplifications imply drastically simpler Boltzmann equations and induce small 
changes to the corresponding decay and production rates. 

Again we solve the set of Boltzmann equations in combination with the initial 
conditions set by preheating and the decay of the gauge degrees of freedom. For our 
key observables we obtain 

Trh ^ 6.1 X 10^ GeV, r/^ ~ 3.7 x 10"^ , r]'^ ~ 3.7 x 10^^ , (69) 
T]'^ ~ 9.7 X 10"^^ , n^h"^ ~ 0.11 . 

With regard to their first two digits, these results for Trh, Vb, Vb ^^h^ are the 
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same as in the full analysis. The result for rf^ is smaller by a factor 2 reflecting the 
missing contribution from the thermal sneutrinos of the first generation. In the upper 
panel of Fig. |5]we present the corresponding comoving number densities. They behave 
very similarly to the original densities in the upper panel of Fig. |2} the only minor 
differences being the following. At early times all densities but the one of the Higgs 
bosons are a bit smaller, at most by a factor of (9(10). In turn, the density of the Higgs 
bosons is technically a bit larger. But the relative change is of O (10~^) and thus not 
visible in Fig. [5j The fact that initially more energy remains in the Higgs bosons has two 
reasons. First, there are now simply less particle species present into which the initial 
vacuum energy could be distributed. Second, particles coupling to the gauge sector are 
produced in smaller numbers after preheating due to the absence of the superpartners 
of the B—L vector boson. A direct consequence of the densities being initially slightly 
smaller is that they become sensitive to the decays of the nonthermal Nf neutrinos 
a bit earlier. The onset of reheating and the inversion of the lepton asymmetry, for 
instance, take place at a^jj ^ 4.2 x 10^ and a ~ 3.2 x 10^, respectively, while these 
events occur later, at a^jj 5.3 x 10^ and a 4.6 x 10^, if supersymmetry is fully 
included. However, as soon as the R and B—L abundances are dominated by the decay 
products of the Nf neutrinos, the differences between the two plots in the upper panels 
of Figs. [2] and [s] begin to vanish. From a ~ 10^ onwards, they are, up to a factor 2 
between the curves for the thermal (s)neutrinos, at or below the percent level. 

It is easy to understand why the omission of the heavy superparticles does not have 



any effect on our final results. According to Eq. (14) the initial energy densities of 
the gauge scalar C, the gaugino A, the higgsino the inflaton as well as the heavy 
sneutrinos Ni are monotonic functions of the Higgs-inflaton coupling A. Setting A to 
its maximal value, A = 10~^, we obtain upper bounds on these densities. 



PA 

Po 



< o (ir^) , ^MMiA 



apH 



Po 



Po 



< O (10-^) . (70) 



apH ^0 apH 

We thus conclude that no matter how the dynamics of the above species look like in 
detail, their influence on the reheating process will always be outweighed sooner or 
later by the decay of the much more abundant Higgs bosons. Ignoring these particles 
does hence not affect the outcome of our calculation. Similarly, we can show that only 
the fermionic decays of the Higgs bosons are relevant for reheating. The ratio of 
sneutrinos to neutrinos increases monotonically with the mass Mi, cf. Eq. (39). 
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Our upper bound on this mass, Mi = 3 x 10^^ GeV, then translates into 

O (10-^) . (71) 



Ni ^ 



The nonthermal sneutrinos can hence also be safely neglected. In conclusion, our 
numerical results in Eqs. ( 70 ) and ( 71 ) substantiate our introductory comment at the 
beginning of this subsection. The essential feature of our scenario of reheating is the 
Higgs boson decay chain, a — )■ Nf — )■ R. From the point of view of the final results for 
the observables, the inclusion of the full supersymmetric particle spectrum is rather a 
matter of theoretical consistency than a numerical necessity. 

Nonsupersymmetric analysis neglecting the gauge multiplet 

Finally, we wish to demonstrate that one is also free to neglect the decay of the gauge 
particles if one is only interested in numerical results for the observables. In addition to 
all massive superparticles we now also exclude the B—L vector boson from our analysis. 
Consequently, particle production in the decay of gauge particles does not take place 
any longer, which simplifies our set of Boltzmann equations once more. This time we 
find for our key observables 

Trh - 6.1 X 10^ GeV, r]B ~ 3.7 x 10"^ , r/^* ~ 3.7 x 10"^ , (72) 
r]%' ~ 9.7 X 10"^^ , fig/i^ ~ 0.11 . 



With regard to their first two digits, these results exactly match those in Eq. (69). The 
lower panel of Fig. [5] displays the corresponding comoving number densities, again to 
be compared with the original densities in the upper panel of Fig. |2} The absence of 
(s) neutrinos of the second and third generation produced through the decay of gauge 
particles now results in a slightly smaller initial lepton asymmetry and, more impor- 
tantly, in drastically shorter N2 3 reheating. While this first stage of reheating still 



lasted until a ~ 11 in our complete analysis, cf. Section 4^, it now comes to an end 
already at a ~ 1.7. Before the onset of A'^i reheating the abundances of radiation, ther- 
mal neutrinos and gravitinos are hence significantly reduced. For instance, at a = 50 
the respective comoving number densities are suppressed by factors of the following 
orders of magnitude, 

B-L : O (10^^) , R, N'^, Nl"^ : O (lO'^) , G : O (lO^^) . (73) 

As before, due to this initial suppression these densities are earlier sensitive to the decay 
of the neutrinos. Now the onset of Ni reheating and the inversion of the lepton 
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asymmetry take place at a\^-^ ~ 1.2 x 10^ and a ~ 2.6 x 10^, which is even earher than 
in our nonsupersymmetric analysis including the gauge multiplet. However, during A^^i 
reheating the differences between the two plots in the upper panel of Fig. [2] and the 
lower panel of Fig. [s] again vanish. From a ~ 10^ onwards, they are at or below the 
percent level. In conclusion, we find that including the gauge degrees of freedom has a 
great impact on the dynamics at early times shortly after preheating, but turns out be 
nonessential when calculating the final numerical results. 



5 Scan of the parameter space 



The value of the Boltzmann equations derived in Section 3.2 is twofold. On the one 
hand, as we have seen in the last section, they are the basis for a detailed time-resolved 
description of the dynamics during reheating. On the other hand, as we will demon- 
strate in this section, solving them in the entire parameter space allows one to study the 
quantitative dependence of our key quantities, Trh, Vb, and VL^h?, on the parameters in 
the Lagrangian. In doing so we will mainly focus on the physical aspects of our results, 
referring the interested reader to Ref . [9] , where we elaborate more comprehensively on 
the technical details of our approach. 

The relevant parameters of our model are the scale of B—L breaking vb-l, the 
heavy neutrino mass Mi, the effective neutrino mass mi, the gravitino mass mg, and 
the gluino mass rrig. Requiring consistency with hybrid inflation and the production 
of cosmic strings fixes the B—L breaking scale, vb-l = 5 x 10^^ GeV, and limits the 



range of possible Mi values, cf. Section 2.4 According to the Froggatt-Nielsen flavour 
model, rhi should be close to fn^ ~ 3 x 10~^ eV. However, in order to account for the 



uncertainties of the flavour model, we vary it between 10~ eV and leV, cf. Eq. (13). 
For the gravitino mass we consider typical values as they arise in scenarios of gauge or 
gravity mediated supersymmetry breaking, 

30 MeV <m^< 700 GeV . (74) 

As for the gluino, we stick without loss of generality to the mass that we used in the 
parameter example discussed in the previous section, rrig = 1 TeV. The generalization 
to different choices for rrig is straightforward, cf. Appendix D in Ref. jo], and simply 
amounts to a rescaling of all values for the gravitino mass. Notice that gravitino masses 
as large as 700 GeV are, in fact, inconsistent with unified gaugino masses at the GUT 
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scale. If the gluino and the bino had the same mass at the GUT scale, the different 
running of the respective renormalization group equations would then entail a mass 
ratio of roughly 6 at low energies. The gravitino which we assume to be the lightest 
superparticle would then have to be lighter than the bino, resulting in an upper bound of 
mg < 170 GeV. We, however, leave open the question whether gaugino mass unification 
takes place at the GUT scale and work in the following with the full gravitino mass 



range specified in Eq. (74). 

At each point of the parameter space defined by the above restrictions we solve 
the Boltzmann equations and record all important numerical results, which we now 



discuss in turn. In Sections |5.1| and |5.2| we study the parameter dependence of the 
reheating temperature and the final baryon asymmetry, respectively. In particular, we 
devote attention to the composition of the asymmetry in terms of a nonthermal and 
a thermal contribution. By imposing the condition that the maximal possible baryon 
asymmetry be larger than the observed one, we identify the region in parameter space 



that is consistent with leptogenesis, cf. the comment below Eq. (64) 



-- Vb + Vb > V°b" - 6.2 X 10"^° 



= rf^ + > r^obs ^ g 2 X 10"^° . (75) 



In Section 5.3 we then turn to the generation of dark matter in the form of gravitinos. 



Requiring the final gravitino abundance to match the observed density of dark matter, 

fig/i^ = figj^/i^ ~ 0.11, (76) 

we are able to derive relations between the neutrino parameters Mi and rhi and the 
superparticle masses mg and rrig. Combining the two conditions in Eqs. (75) and (76), 
we are eventually even able to set a lower bound on in terms of rhi. 



Note that in all plots in this section (Figs. |6| [7] and |8]) the position of the parameter 
point which we investigated in Section |4] is marked by a small white circle. 

5.1 Reheating temperature 

The process of reheating after the B—L phase transition is accompanied by an inter- 
mediate plateau in the decline of the temperature, which determines the characteristic 



temperature scale of reheating. In Section 4.3 we concretized this intuitive notion and 



defined the reheating temperature Trh as the temperature of the thermal bath at the 
moment when the decay of the Nf neutrinos into radiation is about to become efficient. 
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Figure 6: Contour plot of the reheating temperature Trh as a function of the effective neutrino mass fhi 



and the heavy neutrino mass Mi . The reheating temperature is calculated according to Eq. ( 59 ) after 



solving the Boltzmann equations, cf. Appendix [C] for a comparison of our definition of the reheating 
temperature with other common approaches. The thick horizontal gray lines represent the lower and 
the upper bound on Mi , respectively, which arise from requiring consistency with hybrid inflation and 



the production of cosmic strings during the B—L phase transition, cf. Eq. (13 1. The small white circle 
marks the position of the parameter point discussed in Section |4j 



cf. Eq. (59), 



r;^^ (qrh) = H (aRH) , ^RH = T (aRn) • 



In Appendix [C] we argue that this definition is particularly convenient compared to 
alternative approaches because it is not only representative for the temperature plateau 
during reheating, but also associated with a physical feature in the temperature curve. 
Having at hand the solutions of the Boltzmann equations for all allowed values of 



fhi and Mi, Eq. (59) enables us to determine the reheating temperature as a function of 
these two parameters, Trh = Trh (^i, Mi). As the reheating process is solely controlled 
by Higgs and neutrino decays, Trh obviously does not depend on the gravitino or 
gluino mass. In Fig. [6] we present the result of our analysis. We find that, within the 
considered range of neutrino parameters, the reheating temperature varies by almost 
five orders of magnitude. For rhi = 10~^ eV and Mi = 10^ GeV we have, for instance. 
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Trh ~ lO^GeV, while for mi = IQ-^eV and Mi = 10^^ GeV we obtain Trh ~ 3 x 
10^^ GeV. Remarkably, the reheating temperature never exceeds the neutrino mass Mi. 
Instead it is typically smaller than Mi by one or even two orders of magnitude. As the 
ratio Mi/Trh controls the strength of washout process during reheating, we conclude 
that the effect of washout on the generation of the lepton asymmetry is in most cases 



negligible, cf. Section 5^ where we will come back to this observation. 

The reheating temperature increases monotonically with both neutrino parameters, 
fhi and Mi, with the dependence on Mi being much more pronounced than the de- 
pendence on rhi. In the following we will derive a simple semianalytical approximation 
for Trh by means of which this behaviour can be easily understood. A more detailed 
discussion can be found in Appendix C of Ref . |9] . By definition, Trh corresponds to 
the decay temperature of A'^i neutrinos decaying with the effective rate F^^. To first 
approximation, we may thus write 

90 V" c;— - ,„./ 90 



where 7 = 7 (mi. Mi) denotes the average of the relativistic Lorentz factor relating 
to the vacuum decay rate F^^. This first estimate of the reheating temperature 
fails to accurately reproduce our numerical results because of two imprecisions. First, 



Eq. ( 77 ) is based on the assumption that at a = orh the dominant contribution to the 
total energy is contained in radiation. This is, however, never the case. At a = orh 
the decays of the Nf neutrinos have just set in, so that at this time a significant 
fraction of the total energy is hence always still stored in these neutrinos. On top of 
that, for F^ ^ F;^^, which is the case in almost the entire parameter space, the Higgs 
bosons have not decayed yet at a = crh, so that, in the end, they dominate the total 
energy density at the time of reheating. To remedy this first imprecision, we have to 
multiply Eq. (77) by a~^^^, where a = a (mi. Mi) = ptot (o-rh) / Pr («rh)- The second 
imprecision is related to the fact that we do not explicitly solve the Friedmann equation 
to determine the Hubble parameter, but rather calculate it as a/a with the scale factor 



a being constructed as described in Section 3.2 As a consequence of this procedure, H 
does not always exactly fulfill the Friedmann equation. We account for this technical 



imprecision by multiplying Eq. (77) by /3 where 13 = (3 (mi. Mi) relates d/a to the 
exact solution of the Friedmann equation at a = aRH- For appropriate functions a, (3 
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and 7, the reheating temperature Trh can then be written as 

90 



1/2 



r.lxlOiiGeVxa-i/^r'/V'/' 



0.04 eV/ VlOiiGeV 



The dependence of a, /3 and 7 on mi and M\ follows from the solutions of the Boltz- 
mann equations. Restricting ourselves to the region in parameter space in which 
r^yr° > (9(100), we find that /3 and 7 are basically constant. We obtain /3 ~ 0.99 
and 7 ^ 85 with deviations around these values of a few percent. The dependence of 
the correction factor a on mi and M\ is well described by 

mi \ AO^GeV^ 



a ~ 1.2 X lO'^ X — . (79) 

' 0.04eV/ V Ml ' ^ ' 



Such a behaviour directly follows from the interplay of the decay rates and F^. For 
large F^^ and small F^ reheating takes place quite early, at a time when most Higgs 
bosons have not decayed yet. For small F^^ and large Y\ reheating takes place later 
and not as many Higgs bosons are present anymore at a = qrh- The magnitude of 
a is hence controlled by the ratio F^^F^. which scales like rhi/Mi. This explains the 



parameter dependence in Eq. (79). Putting all these results together yields a fitting 
formula for the reheating temperature that reproduces our numerical results with an 
error of less than a percent in almost the entire parameter space, 

1/4 / \ 5/4 



5.2 Baryon asymmetry 

Based on the solutions of the Boltzmann equations we calculate the nonthermal and 



thermal contributions to the final baryon asymmetry, cf. Eq. (63), for all values of the 
neutrino parameters mi and Mi. We present the result of this analysis in Fig. [7j The 
parameter regions in Fig. [t] where the nonthermal and thermal baryon asymmetries 77^* 
and yf^ are consistent with the observational bound r]"^^ are shaded in bright green and 
gray green, respectively. The overlap of these two regions is coloured in dark green. 
In the white patch around mi ~ 0.3 eV and Mi ~ 10^^ GeV the total asymmetry 
Vb = V'b'^V^b is larger than t]'^'^, but neither of its two contributions is. Below the solid 
blue line in Fig. [7] the nonthermal asymmetry dominates over the thermal one. Above 
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Figure 7: Contour plot of the baryon asymmetry rys as a function of the effective neutrino mass fhi 
and the heavy neutrino mass Mi. The baryon asymmetry is calculated according to Eq. (|63]) after 
solving the Boltzmann equations. In the bright green (gray green) region the nonthermal (thermal) 
asymmetry is consistent with the observed asymmetry. In the red region the total asymmetry falls short 
of the observational bound. Below (above) the thin blue line the nonthermal (thermal) asymmetry 
dominates over the thermal (nonthermal) asymmetry. The thick horizontal gray lines represent the 
lower and the upper bound on Mi, respectively, which arise from requiring consistency with hybrid 



inflation and the production of cosmic strings during the B—L phase transition, cf. Eq. (13). The 
small white circle marks the position of the parameter point discussed in Section |4j 



the solid blue line it is the other way around. We conclude that in the part of parameter 
space that we are interested in, the thermal asymmetry is almost always outweighed by 
its nonthermal counterpart. Especially in the region in which leptogenesis is consistent 



with gravitino dark matter, where Mi is typically of O (10^^) GeV, cf. Section 5.3, the 
thermal asymmetry is negligibly small. 

In most of the parameter space the nonthermal asymmetry is insensitive to rhi and 
thus solely controlled by Mi. Only for large values of fhi and Mi it depends on both 
neutrino mass parameters. This behaviour is directly related to the efficiency of the 
washout processes in the respective parameter regions. Let us suppose for a moment 
that washout does not take place. The final nonthermal asymmetry then only depends 
on the total number of (s)neutrinos produced during reheating and the amount of CP 
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violation per (s)neutrino decay. Neither of these two quantities is, however, affected 
by changes in rhi, so that the asymmetry, indeed, ends up being a function of Mi 
only. From this perspective, the insensitivity of r/^* to fhi signals that the effect of 
washout on the generation of the asymmetry is negligible for most values of the neutrino 
parameters. This result is consistent with our findings for the reheating temperature 



and in particular the ratio Mi/Trh as a function of mi and Mi, cf. Section 5.1 To 
see this, note that for temperatures T < Mi the effective washout rate Tyy decreases 
exponentially when raising the ratio Mi/T, 



T < M ■ f - ^1 ^ I 1 ^-Ah/T-pO 

J ^ JWi . i-w - ^J^^ Ni^ j e ' i , 



which readily follows from Eqs. (41) and (42). The fact that Mi/Trh is of (9(10) 
or even larger for most parameter values then explains why the impact of washout 



is typically vanishingly small. In turn, Eq. (81) also illustrates the importance of 
washout at very large values of mi and Mi, for which the ratio Mi/Trh approaches 
values of 0{1). Comparing our results for the reheating temperature and the baryon 
asjTiimetry in Figs. [6] and [7} respectively, we find that washout only plays a significant 
role if Mi/Trh ^ 10 and Mi > 10^^ GeV. Interestingly, the parameter region defined 
by these two conditions covers the entire range of parameters in which the thermal 
asymmetry exceeds the observed asymmetry. 

If washout is negligible, the nonthermal asymmetry can be reproduced to good 
approximation by assuming that all Nf neutrinos decay instantaneously at time ti = 
ts + l/rATj into radiation. The resultant baryon asymmetry is then given by 



nt J^,s i 



(82) 



where denotes the average energy per neutrino. The ratio T /ef^^ is proportional 
to N^_^/Nfi, the number density of Nf neutrinos at the same time when these decay, 
normalized to the radiation number density. It directly follows from the solutions of 
the Boltzmann equations and is well described by 



T 



t=ti 



Together with the expression for ei in Eq. (47) this yields the following fitting formula 
for the nonthermal asymmetry in the case of weak washout, 

/ M \ 
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It reproduces our numerical results for t]'^ within a factor of 2 for most values of Mi. 

The requirement that the maximal possible asymmetry be larger than the observed 
one constrains the allowed range of Mi values. Fig. [7] implies the following lower bound, 

Vb > - 6.2 X 10"^° — > Ml > Mf" ~ 1.7 x 10^° GeV , (85) 

where we have averaged out the slight dependence on rhi. If Mi is chosen below this 
minimal value, the asymmetry falls below the observational bound for two reasons. On 



the one hand, small Mi implies a small CP parameter ei, cf. Eq.(47). On the other 



hand, according to Eq. (83), a small Mi value also entails a small ratio T /efj , i.e. a small 



abundance of (s) neutrinos at the time the asymmetry is generated. The combination of 
both effects then renders the successful generation of the lepton asymmetry impossible. 

The thermal asymmetry has, to first approximation, the same parameter depen- 
dence as the asymmetry generated in standard leptogenesis. It increases monotonically 
with Ml. If Ml is kept fixed at some value Mi > 10^^ GeV, it is largest for rhi values 
of O (10^^) eV. The monotonic behaviour in Mi is a direct consequence of the fact 
that the CP parameter ei scales linearly with Mi. The preference for intermediate 
values of rhi has the same reason as in the standard case. Large fhi corresponds to 
strong washout, at least for the high values of Mi at which the thermal generation 
of the asymmetry carries weight. Small fhi results in a low temperature and a small 
neutrino decay rate such that the thermal production of (s) neutrinos is suppressed. 
Especially in the parameter region in which the thermal asymmetry dominates over the 
nonthermal asymmetry, the expectation from standard leptogenesis ?7§ approximates 
our numerical results reasonably well, 

^^g = |— aphei«:/(mi). (86) 
^ 9*,s 

Here, Kf = (mi) denotes the final efficiency factor. In the strong washout regime, 
•fhi ^ 10^'^ eV, it is inversely proportional to rhi and independent of the initial condi- 
tions at high temperatures p!2] , 

K/(mi) ~ 2 X lO^M — j . (87) 



Combining Eqs. (86) and (87) with the expression for ei in (47), we obtain 



r/f ^ 7.0 X 10-1° 



OAeVVW Ml 



mi J 110^2 GeV 
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In the region in parameter space where ifg > 7]^ this fitting formula reproduces our 
numerical results within a factor of 2. 

Despite these similarities it is, however, important to note that our thermal mech- 
anism for the generation of the lepton asymmetry differs from the standard scenario 
in two important aspects. First, our variant of thermal leptogenesis is accompanied 
by continuous entropy production, while one assumes an adiabatically expanding ther- 
mal bath in the case of standard leptogenesis. Consequently, our thermal asymmetry 
experiences an additional dilution during and after its generation, cf. the comment on 



page |33j Second, our scenario of reheating implies a particular relation between the 
temperature at which leptogenesis takes place, which is basically Trh in our case, and 
the neutrino mass parameters, cf. Section 5.1[ that differs drastically from the corre- 
sponding relation implied by standard leptogenesis. This translates into a different 
parameter dependence of the ratio Mi/T as a function of rhi and Mi, which in turn 
alters the efficiency of washout process and the production of thermal (s)neutrinos from 
the bath in the respective regions of parameter space. In the end, our thermal asym- 
metry therefore rather corresponds to a distorted version of the asymmetry generated 
by standard leptogenesis. As we have remarked above, in the parameter region where 
the thermal asymmetry is larger than the nonthermal asymmetry rf^ hardly deviates 
from 7]^. But as soon as we go to smaller values of rhi and Mi the difference between 
the two asymmetries grows. The minimal value of Mi for which the thermal asymme- 
try is still able to exceed the observational bound, for instance, turns out to be much 
larger in our scenario than in standard leptogenesis. We find an absolute lower bound 
on Ml of roughly 5.1 x 10^^ GeV at an effective neutrino mass rhi — 3.3 x 10^^ eV, 
while standard leptogenesis only constrains Mi to values larger than Mi ~ 10^ GeV. 
Lowering Mi below 5.1 x 10^^ GeV either implies a larger ratio Mi/Trh or a larger 
effective neutrino mass mi, cf. Fig. |6] In either case the thermal asymmetry is reduced 
so that it drops below the observed value. 

In conclusion, we emphasize that the generation of the lepton asymmetry is typ- 
ically dominated by the decay of the nonthermal (s)neutrinos. Only in the param- 
eter region of strong washout, which is characterized by a small ratio Mi/Trh, the 
nonthermal asymmetry is suppressed and the thermal asymmetry has the chance to 
dominate. Related to that, we find that the viable region in parameter space governed 
by the nonthermal mechanism is significantly larger than the corresponding region 
for the thermal mechanism. Independently of mi, the neutrino mass Mi can be as 
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small as M™™ ~ 1.7 x 10^" GeV, which is an order of magnitude below the bound of 
5.1 X 10^^ GeV, which one obtains in the purely thermal case. 

5.3 Gravitino dark matter 

The final abundance of gravitinos Q^h"^ depends on three parameters: the reheating 
temperature Trh as well as the two superparticle masses mg and rUg. A key result 
of our reheating scenario is that Trh is determined by the neutrino mass parameters 
fhi and Mi. As we keep the gluino mass fixed at ITeV, the gravitino abundance thus 
ends up being a function of fhi, Mi and mg. Based on the solutions of the Boltzmann 



equations we calculate Q^^h"^ according to Eq. (67) for all values of these three masses. 



By imposing the condition that gravitinos be the constituents of dark matter we can 
then eliminate one of the free mass parameters, for instance the neutrino mass Mi, 

n^h^ (mi. Ml, mg) = ^dm^^ — ^ Mi = Mi (mi, mg) . (89) 

The physical picture behind this step is the following. For given mg, the reheating 
temperature has to have one specific value so that the abundance of gravitinos comes 
out right. Each choice for mi then implies one particular value of Mi for which this 



desired reheating temperature is obtained. Solving Eq. (89) for Mi yields this value as a 
function of mi and mg. The corresponding reheating temperature follows immediately, 

Trh = Trh (mi. Ml (mi, m^)) — )■ Trh = Trh (mi, mg) . (90) 

In summary, combining the requirement that gravitinos make up the dark matter with 
the fact that the reheating temperature is determined by neutrino parameters allows 
us to infer relations between these neutrino parameters and superparticle masses. The 



lower bound on Mi induced by leptogenesis, cf. Eq. (85), can then be translated into a 
constraint on the mass parameters mi and mg. 

Vb = Vb {mi, Mi (mi, mg)) > 7]°j^' — > Vb = Vb {mi, mg) > r/^*'' . (91) 

We present our results for the functions Mi (mi,mg) and Trh ("^i,"^^) in the two 
panels of Fig. |8[ respectively. Furthermore, we indicate in both plots the constraint 
arising from the requirement of successful leptogenesis. 

We observe the following trends in the two plots of Fig. |8} Both quantities. Mi and 
Trh, show a stronger dependence on the gravitino mass than on the effective neutrino 
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Figure 8: Contour plots of the heavy neutrino mass Mi (upper panel) and the reheating temperature 
Trh (lower panel) as functions of the effective neutrino mass mi and the gravitino mass mg such 



that the relic density of dark matter is accounted for by gravitinos, cf. Eqs. (89) and (90). In the red 



region the lepton asymmetry generated by leptogenesis is smaller than the observed one, providing us 
with a lower bound on the gravitino mass in dependence on fhi. The colour code is the same as in 
Figs. |6] and [7| The small white circle marks the position of the parameter point discussed in Section |4] 
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mass. For fhi < 10^^ eV the reheating temperature is almost completely insensitive to 
fhi. The neutrino mass Mi slightly increases when lowering the value of fhi. For large 
values of the effective neutrino mass, fhi > 10~^ eV, exactly the opposite is the case. 
Ml does not depend on fhi anymore and Trh slightly rises when increasing fhi. In the 
following we will construct semianalytical approximations for Mi and Trh which will 
allow us to get some intuition for this behaviour. The final gravitinos abundance Q^h"^ 
can be parametrized in the following way, cf. Appendix D of Ref. (ol for details. 



eCi 



RH 



IQio GeV 



Co 



100 Gev) 



100 GeV 



ITeV/ 



(92) 



Here, the two coefficient functions Ci,2 = C'1,2 (^rh) subsume all factors contributing 
to ^?^h? that can be taken care of analytically. 



Ci = 10^* GeV 



90 



1/2 



ISgt (Trh) 
g^g^ (/xo) Mp 



log 



-^RH 



^1 (Trh) 



+ 0.8846 



C2 



if^o) 



(93) 



WOgt (Trh) ' 

They both depend only very weakly on the reheating temperature, so that for our 
purposes it will suffice to treat them as constants, Ci ~ 0.26 and C2 — 0.13. The factor 
e parametrizes all effects that cannot be accounted for analytically in the derivation of 
Eq. (92), i.e. the amount of energy in radiation at a = Orh, the ratio ^q/H at a = orh 
as well as the increase in the comoving number densities of gravitinos and radiation 
after a = orh- In principle it depends on all mass parameters, in practice after solving 
the Boltzmann equations we find that it is mainly controlled by fhi, 

lO-^^eV^ " 



where the exponent c is c 



e (mi) ~ 1.2 

mi 

0.21 for fhi > 10~^eV and c 



(94) 



-0.01 for fhi < 10-3 eV. 
We insert our results for Ci^2 and e into Eq. (92), set Q^h'^ to Q^^h'^ and solve for Trh, 
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-1 



(95) 



The corresponding expression for Mi can then be obtained by exploiting Eq. (80), 

-4/5 



Ml ~ 7.2 X 10^" GeV 
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0.13 



100 Gev) 



100 GeV \ 



(96) 
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where the exponent d is given as 4c/5— 1/5 so that d ^ —0.03 for rhi > 10"'^ eV and 
d ~ —0.20 for rhi < 10~^eV. These two fitting formulae reproduce our numerical 
results with deviations of (9(10%) and nicely illustrate the different dependence of Trh 
and Ml on fhi for small and large values of mi, respectively. As expected, they show 
that the dependence on fhi is always very mild and solely stems from the factor e, 
i.e. corrections beyond the purely analytical result for Vt^h"^. If we were to omit these 
corrections and set e to 1, the reheating temperature required for gravitino dark matter 
would be a function of only, Trh = ^rh (^g), in accordance with the fact that 
the only parameters entering the gravitino production rate are the masses of the 
gravitino and the gluino. 

The relation between the gravitino mass and the neutrino parameters rhi and Mi 
translates the lower bound on Mi imposed by the requirement of successful leptogenesis, 
cf. Eq. (85), into a lower bound on mg. As we can read off from Fig. [sj mg must be 
at least of C(10) GeV to obtain consistency between leptogenesis and gravitino dark 
matter. In fact, the bound on slightly varies with fhi. For rhi values between 
10~^eV and 10~^eV it monotonically increases from roughly 7 GeV to 17 GeV, from 
fhi ~ lO^^eV onwards it remains at mg ^ 17 GeV. For such low gravitino masses the 
first term in the brackets on the right-hand side of Eq. (96) is negligibly small so that 
the fitting formula for Mi can be easily solved for mg, 

o -;cc"( X"( V''"° (97) 



Imposing the condition that Mi be larger than M™™ ~ 1.7 x 10^° GeV, cf. Eq. (85), 
provides us with an analytical expression for the lower bound on mg, 

1/4-c 



mg>m|-^16Gev(^^) \ (9J 



This estimate reproduces our numerical results with a precision at the level of O (10%). 
Physically, the connection between the bounds on mg and Mi is the following. For 
gravitino masses below O (10) GeV, a reheating temperature Trh ^ O (lO*^"^) GeV is 
required to avoid overproduction of gravitinos. According to our reheating mechanism 



^In physical terms this means that for small gravitino masses mainly the goldstino degrees of 
freedom of the gravitino rather than its transverse degrees of freedom are excited. Cf. Ref. [9] for a 
detailed discussion on how the interplay between the two different production modes of the gravitino 
is reflected in our results for Mi and Trh- 
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such low reheating temperatures are associated with comparatively small values of the 
neutrino mass, Mi < O (10^°) GeV. The low temperature and low mass then entail a 
small abundance of (s) neutrinos at the time the asymmetry is generated and a small 



CP parameter ei, cf. Eqs. (83) and (47), respectively. Both effects combine and result 
in an insufficient lepton asymmetry, rendering dark matter made of gravitinos with a 
mass below O (10) GeV inconsistent with leptogenesis. 

In conclusion, we find that our scenario of reheating can be easily realized in a 
large fraction of parameter space. The two conditions of successful leptogenesis and 
gravitino dark matter, in combination with constraints from hybrid inflation, allow us to 
interconnect parameters of the neutrino and supergravity sector. In particular, we are 
able to determine the neutrino mass Mi and the reheating temperature Trh as functions 
of the the effective neutrino mass mi and the gravitino mass mg. Furthermore, the 
consistency between all ingredients of our scenario indicates preferences for Mi and 
Trh; namely Mi values close to 10^^ GeV and Trh values close to 3 x 10^ GeV. Finally, 
we obtain a lower bound on the gravitino mass of roughly 10 GeV. 



6 Conclusion and outlook 

A phase of false vacuum of unbroken B—L symmetry at the GUT scale can account for 
the observed acoustic peaks in the cosmic microwave background via hybrid inflation. 
Subsequent tachyonic preheating, followed by the decay of heavy gauge and Higgs 
particles and heavy neutrinos sets the initial conditions of the hot early universe. We 
have studied the B—L breaking phase transition for the full supersymmetric Abelian 
Higgs model and given a detailed time-resolved description of the reheating process 
taking all (super)particles into account. The competition of cosmic expansion and 
entropy production leads to an intermediate plateau of constant temperature, during 
which baryon asymmetry and gravitino dark matter are produced. 

The initial conditions of the thermal phase of the universe are determined by the 
parameters of the fundamental Lagrangian, i.e. the masses and couplings of elementary 
particles. Likewise, the constant plateau temperature is fixed by neutrino parameters. 
The temperature scale of reheating is hence no longer an unknown cosmological param- 
eter, but rather an effective quantity that is determined by mass parameters that can in 
principle be measured in experiments. The consistency of hybrid inflation, leptogenesis 
and gravitino dark matter restricts the parameter space. For a gluino mass of 1 TeV 
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we find a lower bound on the gravitino mass of about 10 GeV. Tlie order of magnitude 
of Ml, the mass of the lightest of the heavy neutrinos, is 10^^ GeV. For a wide range 
of light neutrino masses this results in a reheating temperature of order lO^--^° GeV. 

We point out that lowering the scale of B—L breaking would significantly weaken the 
bound on the gravitino mass. If B—L breaking is unrelated to hybrid inflation and takes 
place at a scale vb-l ~ 10^^ GeV, the gravitino could have a mass of O (100) MeV |9j. 
Similarly, for a lower B—L scale reheating would occur at a higher temperature because 
of faster Higgs decays. This would result in a stronger washout of the lepton asymmetry 
generated in (s)neutrino decays. Small vb-l hence implies an upper bound on the 
effective neutrino mass fhi of about 0.1 eV [9j. In this paper we have demonstrated that, 
if B—L is broken at the GUT scale, this restriction does no longer apply, rendering the 
proposed reheating mechanism viable for all reasonable masses of the light neutrinos. 

Tachyonic preheating is a complicated nonequilibrium process, which requires fur- 
ther theoretical investigations. A remarkable result of this work is that the final baryon 
asymmetry and dark matter density are rather insensitive to many of the related the- 
oretical uncertainties, such as the details of the production and relaxation of cosmic 
strings. For instance, even if 50% of the false vacuum energy density is initially stored 
in strings, they quickly loose most of their energy and the effect on the final baryon 
asymmetry and dark matter abundance is negligible. This robustness is due to the fact 
that after all most of the vacuum energy density is transferred to heavy Higgs bosons 
whose slow decays, via heavy neutrinos, dominate the reheating process. 

Throughout our analysis we have assumed that the gravitino is the lightest super- 
particle. However, the proposed mechanism for the ignition of the hot early universe 
also works if the gravitino is very heavy with a neutralino as LSP. In this case ordinary 
WIMP dark matter can be nonthermally produced from gravitino decays. Consistency 
of hybrid inflation, leptogenesis and dark matter density then leads to constraints on 



gravitino and neutralino masses. In Ref. 53 we give a detailed description of this 
alternative scenario. Further important questions concern the effect of the infiaton on 



tachyonic preheating 54 , 55 and possible modifications of superpotential and Kahler 
potential of the symmetry breaking sector in connection with the detailed description 
of the cosmic microwave background, which will be discussed elsewhere. 
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A The supersymmetric Abelian Higgs model 
in unitary gauge 

In this section, we present the full supersymmetric Lagrangian of the Abelian Higgs 



model in unitary gauge, following the notation of Ref. 56 . Our starting point is an 
arbitrary superpotential W, given in terms of chiral fields whose scalar and fermionic 
components are denoted by (pi and tpi, and the canonical Kahler potential 

ir = $tgP.v^.^ p^ = 2gq,, (99) 

where g is the gauge coupling and is the U{1) gauge charge of ^j. V denotes the 
f/(l) vector superfield, 

V = C + iex-iOx+ ^Oe (M + iN) - U)e (M - iN) - ea^'OA^ 

+ ^oee (^^ + l^T^'d.x^ - loee (^i + \o^'d,t^ + hem {d + ]pc^ , (loo) 

containing the scalar degree of freedom C, the fermionic components x ^"^^ the vector 
Afj_ as well as the auxiliary fields D, M and N . In the Wess-Zumino gauge, which we 
will not use, one has C = 0, x = and M = A = 0. 

The supersymmetric Lagrangian can be derived in the standard manner by cal- 
culating D- and F-terms of Kahler potential and superpotential and eliminating all 
auxiliary fields. In order to obtain fields with canonical mass dimension we perform 
the rescalings 

^c^c, -^x-^x, (101) 

where p corresponds to one specific, conveniently chosen pi and v is an arbitrary nonva- 
nishing mass scale. In the following, we will promote f to a time-dependent function. 
Note, however, that our discussion also applies to the even more general case of a fully 
spacetime-dependent scalar field v = v {t^x). 
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After some calculations, including several integrations by part, one finds the La- 
grangian. 



r _ rkin I /'gauge , Mcrm t/ \l \ C 

~ *-WZ "T ''-WZ "T ^WZ — yp — Vd -r ^non-WZ ; 



(102) 
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(103) 



Evaluating the exponential functions in >Cwz to leading order in pi V2C/{pv) yields the 
familiar Lagrangian in Wess-Zumino gauge; the remaining terms, collected in >Cnon-wz, 
represent additional terms involving the gauge degrees of freedom C and x- 
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The symmetry breaking sector defined in Section |2| cf. Eq. ([T]), contains the super- 
fields 5*1, 5*2 and $ with B—L charges qs^ = —Qsi = ^5 = 2 and g$ = 0. In unitary 
gauge, cf. Eq. ([s]), one has 

Si,2 = ^S' . (104) 
With p = ps = 2gqs, S' = (s, s), and $ = (0, 0), one now obtain 



/'kin 



cosh {d^s*d^s + iSa^'d^.s) 



cosh 



= sinh 




^ (is*d^'s - isd^s* + Id^s) 



L 2 



^ I -^s*sf H — "\/A</)ss + \/As0s + h.c. 



o 



cosh 



(105) 



(106) 
(107) 
(108) 
(109) 



and 
-^non-wz = sinh 



.2^2 f 2f2 



s*xcr^(9„s + h.c. 



+ cosh 



+ 



y L ^ ^ V 2t; 



;iio) 



72: 



+ + h.c. 



+ h.c. 



The ground state of the theory corresponds to |sp = Identifying the mass 

scale V with the time-dependent vacuum expectation value of the Higgs field in the 
''For notational convenience, we have omitted the prime on the complex scalar Higgs boson s. 
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broken phase, v = v{t) = -^{a''^{t,x))]l'^ , which approaches vb-l at large times, the 
Lagrangian £non-wz yields kinetic terms for C and x aiid a mass term for x aiid ^. 



The mass terms for and C are contained in Eqs. (106) and (109), respectively. As 



expected, in unitary gauge the vector field V describes a massive vector multiplet 56 
Shifting s around its expectation value, s ^ v{t) + + it), one reads off the masses 
given in Eq. ([6]). 

Note that due to the time-dependence of f, the kinetic term for C in Eq. (110) 
yields a contribution to the mass nic- In the main part of this paper, we omit this 
term for two reasons. First, it is much smaller than the contribution to mc obtained 



from Eq. (109) throughout the preheating process and hence the latter governs the 



production during tachyonic preheating. Second, as we show in Section 4^, our final 
results prove insensitive to the dynamics of the gauge sector and we can hence ignore 
this technically rather complicated contribution. 



B CP violation in 2 — > 2 scattering processes 

To calculate the lepton asymmetry consistently to first order in the CP violation pa- 
rameter e, 2 — )■ 2 scattering processes involving an (anti-)(s)lepton in the initial and 
final state must be considered. Scatterings with an on-shell neutrino in the interme- 
diate state are already included in decay and inverse decay processes. We are hence 
left with the task to calculate the off-shell contribution of these processes. For the 



nonsupersymmetric case, this was discussed in Refs. |49j and 57 . Here we explain the 
supersymmetric case. We first study the CP-violating contribution of the full 2 — 2 
scattering processes and will see that this vanishes to 0{{h'^Y). Hence to this or- 
der in the Yukawa coupling, the CP-violating off-shell contributions can be added by 
subtracting the corresponding on-shell contributions. 

The right-hand side of the integrated Boltzmann equation is given by the interaction 



density 7 = gxi^n) ^ f d^pCx, cf. Eqs (20) and (23). For distinct final and initial 



states, this is related to the corresponding S'-matrix elements 

J2^{^^f) = J2\SFI\'fI, (111) 

where the summation over the lower case letters on the left-hand side runs over dif- 
ferent particle species and the summation over capital letters on the right-hand side 
additionally includes the summation over all internal degrees of freedom as well as 
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phase space integrals for all initial and final state particles. Considering the case of 
2 — 7- 2 scatterings in the Boltzmann equation for the lepton asymmetry, the initial and 
final states of interest are {i,f} G {iH,iH,iH,iH}. The internal degrees of freedom 
are helicity, weak isospin and flavour. // denotes the phase space distribution function 
of particle species i. 

Using this notation, we now consider the CP-violating contributions of the full 
2 — )■ 2 scattering processes, 

E ^ /) - ^(^ ^ /)] = E O^i^/i V/ - \SFi\'fi] 

i,f I,F 

= [l^Fil'^ + I'S'f/P - \Sfi\'^ - \Sfi\'^] fi (112) 

I,F 

= J2i^-i]fi + oi{hr) = oi{hr). 

I 

The bar indicates CP conjugation and // = // are the phase space distributions of 
the light MSSM (anti-)particles in thermal equilibrium. Here in the second line of 



Eq. ( 112 ), we extended the summation over the final states to include the lepton number 
conserving processes. These can be grouped in pairs of CPT conjugates and hence, due 
to CPT invariance, yield a vanishing contribution in total. In the third line, we exploit 
the unitarity of the 5* matrix, i.e. that the summation over all possible final states 



yields 1. Since however in Eq. (112) the sum runs only over all possible two-particle 
final states, we obtain corrections caused by neglecting multi-particle final states. For 
off-shell intermediate states these corrections are of 0{{h'^)^) [57], however close to the 
resonance pole they are enhanced to C((/i^)^) 

Concluding, we find that the CP- violating contributions of the 2—7-2 scattering 
processes involved in the production of the lepton asymmetry vanish, with corrections 
of 0{{h'')'^). Hence the on- and off-shell contributions cancel each other and we can use 
the usual 'recipe' of replacing the CP-violating contributions of the off-shell (s)neutrino 
decays by the negative of the respective on-shell contributions, i.e. 

E E ^(^^' ^ /) = - E E ^(^"'^ ^ /) + ^((^^)') ' (113) 

/a fa 

where a is a flavour index. Note that looking at this line of argument closely, this 
argument holds separately for neutrinos and sneutrinos because of distinct sets of initial 
and final states, but the summation over flavour and lepton/slepton is unavoidable. 
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C Definition of the reheating temperature 



Apart from the definition of the reheating temperature employed in this work, cf. 



Eq. (59), there are alternative ways to define the reheating temperature. For instance, 
we could use the temperature at the beginning (a = aj^jj) or the end of reheating (a = 
^rh) '^^ temperature when half of the total available energy has been transferred 
to radiation (a ~ as for the parameter example discussed in Section |4]). In either case, 
although the respective value for orh may significantly vary, thanks to the temperature 
plateau during reheating the resulting reheating temperature would not change much. 
For the parameter point investigated in Section [4], we find 

^i^^l.5, ZKh1^^_ (,,4) 

^RH ^RH 2.5 Trh 3.0 

Our definition of the reheating temperature may hence be regarded as a compromise 
between several more extreme approaches. But more important than that, it picks 
up on a physical feature that other definitions would miss. In Fig. |4] we observe that 
the temperature declines less during the first part of reheating, a^jj < a < orh, than 
during the second part, orh < « < Orh- '^^^ stage of Ni reheating evidently splits 
up into two phases, during the first of which the temperature is basically constant, 
whereas during the second one the temperature slightly decreases. The reason for this 
substructure in the temperature plateau is the following. As soon as the Nf neutrinos 
decay more efficiently their comoving number density starts to grow slower than a^/^. 



This diminishes the production rate of radiation. According to Eq. (58), a constant 
temperature can then no longer be maintained. The advantage of our definition for 
Trh now is that we read it off the curve in Fig. |4] at exactly that value of the scale 
factor at which the transition between these two phases of A^^i reheating takes place. 
Our definition thus yields a temperature that is both representative as it mediates 
between several more extreme values and especially singled out as it is associated with 
a prominent feature in the temperature curve. 

For completeness, we should however mention that for other parameter choices 
this picture may change. If the Higgs decay rate F^ is, for instance, larger than the 
neutrino decay rate F;^^, which can for example be achieved by going to lower values 
of the B—L scale, the scaling behaviour of the Nf number density changes when the 
neutrino production efficiency begins to cease and not when the decays of the neutrinos 
themselves set in. The slight kink in the temperature plateau is then located at a ^ 
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which is in this case before the decay of the Nf has become fully efficient. But the 



definition of the reheating temperature in Eq. (59) remains reasonable nonetheless. 
After all, if > T^^, the bulk of the total energy is first almost entirely accumulated 
in neutrinos before it is passed on to radiation. The energy in radiation thus receives 
its major contribution just when these neutrinos decay with a sufficient efficiency. The 
characteristic temperature at the time when this happens is then again obtained from 



Eq. (59). Further details on the reheating temperature in regions in parameter space 



in which > Vf^^ can be found in Ref. |9 
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